Tailoring Three-Point Functions and Integrability 



Jorge Escobedo"''', Nikolay Gromov"^, Amit Sever" and Pedro Vieira" 



O 



^ Department of Physics and Astronomy & Guelph- Waterloo Physics Institute, 
University of Waterloo, Waterloo, Ontario N2L 3G1, Canada 



Perimeter Institute for Theoretical Physics 
Waterloo, Ontario N2L 2Y5, Canada 
j escob , amit . sever , pedrogvieiraOgmail . com 

Oh 
< 

o 

^ " King's College London, Department of Mathematics WC2R 2LS, UK & 

St. Petersburg INP, St. Petersburg, Russia 
nikgromovQgmail . com 

Abstract 

> 

in 

^ We use Integrability techniques to compute structure constants in A/" = 4 SYM to leading 

^ order. Three closed spin chains, which represent the single trace gauge-invariant operators in 

^ J\f — 4: SYM, are cut into six open chains which are then sewed back together into some nice 

O pants, the three-point function. The algebraic and coordinate Bethe ansatz tools necessary 
for this task are reviewed. Finally, we discuss the classical limit of our results, anticipating 

. ^ some predictions for quasi-classical string correlators in terms of algebraic curves. 
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1 Introduction 



Solving interacting conformal field theories in Ad with a large expansion will have a deep 
impact in our understanding of Nature. M = A super Yang-Mills (SYM) seems to be, 
excitingly, the harmonic oscillator of gauge theories in four dimensions. Its full solution will 
dramatically improve our understanding of particle theories such as QCD. 

The fundamental objects in conformal field theories are two- and three-point functions 
of local gauge-invariant operators. Knowing them we can in principle construct any higher- 
point function by gluing these building blocks together. Two-point functions are greatly 
understood in A/" = 4 SYM, largely due to the existence of integrability [1, 2]. In this paper 
we shall focus on the study of planar three-point functions, or structure constants, using the 
underlying exactly solvable structures of these theories. We will illustrate our methods at 
weak coupling in A/" = 4 SYM. 

For some interesting and inspirational works on three-point function in TV = 4 SYM at 
weak coupling see [3, 4, 5]. In particular, [3] introduces the physical picture of cutting and 
gluing spin chains which we elaborate on below and [4] emphasizes the usefulness of the 
algebraic Bethe ansatz techniques for computing scalar products of quantum spin chains, 
which turns out to be very relevant for this problem. 

In any conformal field theory (CFT), one can choose a basis of local operators such that 
their two-point functions are given by 

mx)oM)=^^T^ . (1) 

I I ^ 

where Aj are their conformal dimensions and Mi are normalization constants that may be 
set to one. The correlation function of three such local operators is restricted by conformal 
symmetry to be of the form 



{0.:(x<)0,(-,)0.(x.)> = (2) 

l-^ljl l-^jfel l-^Kll 

where Cijk are the structure constants. In what follows we will compute these structure 
constants in the planar limit of A/" = 4 SYM to leading order in l/N^ and A, i.e. the number 
of colors and the 't Hooft coupling. The first nonzero structure constants in the 1/Nc 
expansion scale as 1/Nc. These arise in the three-point function of single trace operators. 
Therefore, our operators will all be of the form Ti{ABC . . . ) where A, B, C are A/" = 4 fields. 
These are the most basic correlation functions building blocks.^ 

The corresponding single trace structure constants have a perturbative expansion of the 
form 

NcC,^^ = c^ + \^ + >?<^ + ... (3) 

^One would expect that in order to Bootstrap the four-point function of single trace operators, the 
knowledge of the three-point functions between two single trace operators and one higher trace operator 
would also be required. Remarkably, this does not seem to be the case; the four-point function seems to be 
re-constructable from the three-point function of single trace operators alone [6] . 
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and we will only consider the first term in this expansion c-jl- At A = many single trace 
operators have the same dimension, which is roughly the number of fundamental fields in 
the operator. That huge degeneracy is lifted at one loop. Therefore, to correctly identify the 
coefficient cf^l in the expansion (3), we have to use those linear combinations of single trace 

operators that have definite one-loop anomalous dimension Aj = A^-"'' + X^y^^^-'^ 

We are therefore lead to the following picture (see figure la). The three single trace 
operators of lengths Li, L2 and L3 are contracted by free propagators. Since each propagator 
connects two fields, Li + L2 + L3 must be an even number. The number of free contractions 
between Oi and Oj is (Lj + Lj — Lk)/2. These propagators automatically reproduce the 
factor l/\xij\^^'^^^~^'' in (2), with Aj = the free dimension A-^''. In the planar limit, these 
propagators are all color neighbors. The tree-level structure constant is then given by the 
sum over all such contractions, normalized by the two-point functions. This paper is devoted 
to the study of this interesting combinatorial problem. 

The three operators have definite one-loop anomalous dimension and are given by a linear 
combinations of single trace operators. Linear combinations of single trace operators can be 
represented as states |\E'j) on a closed spin chain. For example, the member of the Konishi 
multiplet composed from two complex scalars Z and X can be represented in a spin half 
chain by the state 

K oc Tt{ZZXX) - Tt{ZXZX) = |ttU) - Ititi) • (4) 

In the spin chain language, the structure constant is constructed by going through the 
following steps (see also figure lb) 

1. Starting with three closed chains in the states l^&i), |^fc) and choosing a cyclic 
ordering (z, j, k), 

2. Breaking the ith closed chain into left and right open subchains of lengths {Li + Lj — 
Lk)/2, {Li + Lk — Lj)/2 and doing the same for the other two closed chains, 

3. Expressing the closed chain state as an entangled state in the tensor product of the 
two subchains Hilbert spaces with the lengths indicated in the previous point. I.e., 

4. We want to Wick contract the operator corresponding to the state \'^i)r with the 
operator corresponding to the state |\E'j+i)i. The Wick contraction is obtained from 
the spin chain contraction of a ket and bra states ^(^il^I'i+i)/ after a flipping operation 
|\E'i);® |\E'j)r — )■ |\I'j)/® r(^j|, see figure lb. This fiipping operation maps the ket states 
in the right subchain into bra states with 1) reversed spin chain sites, 2) same wave 
function (not conjugated), 3) same charges ^, 

5. Normalizing the three external states. 

■^This is nothing but the standard textbook degenerate perturbation theory in QM and needs to be taken 
into account. See e.g. [7] and [8] for a discussion of the importance of this point. 

•^This point might be a bit confusing at first. It is explained in greater detail in section 3 and was 
previously considered in [3]. To help the curious reader let us give an example of this map of ket into 
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(a) 



(b) 



Figure 1: (a) The planar tree- level contraction of three single trace operators in the double line 
notation. The diagram has a pair of pants topology (sphere with three punctures), (b) In the spin 
chain picture, each of the three single trace operators corresponds to a state on a closed chain. The 
closed chains are cut into right and left open chains where the external states are represented. The 
three states are sewed together into the three-point function by overlapping the wave functions on 
each right chain with the wave function on the left subchain of the next operator. 



The resulting structure constant can then be obtained from brute force contractions of 
the states and is given by 



JO) 
-'ijk 



LiL.Lk E r(^fcj^ij« r(^,J^'fc.)i 
a,b,c 



(5) 



Here, the factors of Li, Lj and arise from summing over all the ways of cutting open the 
closed chains, before gluing them together. 

What we have said so far applies to the planar tree-level contraction of any single trace 
operators, whether they have a definite one-loop anomalous dimension or not. What we have 
to do, however, is to consider the operators with definite anomalous dimension. Generically 
these operators are some linear combinations of a huge number of single-trace operators. 
Even at tree-level the direct calculation of these contractions is a very complicated problem, 
especially for long operators. One may hope, however, that the integrability hidden in the 
dilatation operator may help to simplify this problem. Indeed, the one-loop anomalous di- 
mension matrix is represented by an integrable spin chain Hamiltonian [2] . We can therefore 
use integrability techniques to compute the normalization of the external eigenstates, called 
Bethe eigenstates. Moreover, at one loop the spin chain Hamiltonian only acts locally on the 
chain. As a result, when we decompose the external state as an entangled state in the two 

hm-.jZZXZX) ^ {XZXZZ\ = {\XZXZZ)y 7^ {\ZZXZX))''. Note also that in our notations {Z\Z) = 
{Z\Z) = 1, {Z\Z) = and the spins in the kets and bras are ordered from left to right, i.e. {ZX\ZX) = 1 
while {ZX\XZ) = 0. 
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subchains, each of the subchain states still has the same local form as an eigenstate. A state 
of that local form is called a Bethe state. Therefore, even though the subchain states are 
not eigenstates, we can still use integrability techniques to compute their overlaps. We will 
show that the integrability techniques applied to this problem are much more efficient than 
a brute force calculation and in particular allow us to make computations for asymptotically 
long operators, which is otherwise impossible. 

To summarize, the tools we need in order to compute the tree-level structure constants 

are 

• The spin chain Bethe eigenstates. 

• The decomposition of an external eigenstate into an entangled state on the direct 
product of the two subchains. We shall denote this decomposition procedure by cutting. 

• Once the state is cut into two we need to flip one of its halves from a ket into a bra, 
see figure lb. 

• The overlaps and norms of Bethe states. We denote the overlapping computations by 
sewing. 

The necessary tools for the cutting, fiipping and sewing of eigenstates into a three-point 
function are explained in section 3. For completeness, the Bethe Ansatz in JV = 4 SYM and 
the construction of the corresponding wave functions is reviewed in section 2. The three- 
point functions are built using the three steps mentioned above in section 4. In this section 
we also discuss some interesting limits such as the BMN limit or the classical limit. Section 
5 contains the conclusions and open problems. Appendices A,B,C,D contain supplementary 
details while appendix E contains examples of structure constants at one loop order and is 
intended to be self-contained. 



2 Bethe Ansatz. Review and Notation 

In this section we will review some basic facts about two-point functions in A/" = 4 SYM, the 
map between operators and quantum spin chain states of an integrable spin chain and the 
coordinate and algebraic Bethe ansatz which can be used to efficiently study the spectrum 
problem at one loop. One of the goals of this section is to set the notation for the following 
sections. The expert reader is certainly familiar with the content of this section and can 
safely skip it. 

In this paper we will mostly consider operators made out of two complex scalars. One 
such example is the Konishi operator presented in (4) made out of the scalars Z and X. As 
mentioned above, we can represent such operators as states in an SU (2) spin chain. One of 
the scalars is thought of as being a spin up - the Z field in the example (4) - while the other 
scalar is thought of as being a spin down - the X field in the example (4). More generally we 
represent the operator made of L scalar fields Z by a ferromagnetic vacuum state of L spins 
up [9] . Operators with N scalar fields X and L — N fields Z are represented by fiipping N 
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of those spins. These spin flip excitations are called magnons. In the example (4) we have 
L = 4 and = 2. A generic state will be of the form 

1^) = ^ tp{ni, . . . ,nN)\ni, . . . ,nN) (6) 

l<ni<n2<---<njv<-L 

where the ket \ni, . . . ,71^) stands for the state with spins down at positions ni,n2, . . .un^ 
while ip{ni, . . . , n^v) is the wave function, which we will flx below. In what follows, whenever 
we will refer to a wave function, we will mean the wave function in this local base. 

The operators with deflnite anomalous dimensions are the eigenvectors of the one-loop 
mixing matrix H (by deflnition) . The anomalous dimensions are the corresponding eigenval- 
ues (again by deflnition). This mixing matrix can be computed in perturbation theory. In 
the spin chain language, the mixing matrix is represented by a local spin chain Hamiltonian 
of the form [2] 

A ^ 

H=—J2iKn+l-^n,n+l). (7) 
n=l 

The permutation operator P„,„+i acts on the spins at positions n and n + 1 swapping them. 
The identity In,n+i does nothing and L + 1 = 1. The energy spectrum of this Hamilto- 
nian gives us the anomalous dimensions of the operators. We encourage the reader who is 
not familiar with this language to check that the state (4) is indeed an eigenvector of this 
Hamiltonian with anomalous dimension 'Jk = 3A/47r^. 

We will now review how to construct the wave functions iplni, . . . , njy) which diagonalize 
the spin chain Hamiltonian (7). For pedagogical references on the coordinate Bethe ansatz, 
see [10, 11, 12]. 



2.1 Coordinate Bethe ansatz 



For a single magnon we diagonalize the translation invariant Hamiltonian by going to Fourier, 
il^{ni) = e*^!"^^ The energy of this state, also called the dispersion relation, is given by 
7 = with 

e(p) = ttV ^^^^ 



— n sm - 
27r2 2 



For two particles we write ip{ni,n2) = e^^'^"'^'^^^^"^ + S {p2, Pi)e^^^^'^~^'^^^"'^ . The relative coef- 
flcient between the two plane waves is the amplitude for incoming momenta {pi,P2} to be 
exchanged into {p2,Pi}- In other words, it is the two body S-matrix. By explicitly acting 
with H on this state we can read the energy of the state 7 = e(pi) + e(p2) and compute the 
S-matrix, 

1 Pot 2 - 1 rot 2^ 7 

s{p,P') = i I ; • (9) 

I cot I - I cot ^ - 2 

Now, in a 1+1 dimensional elastic scattering process between two identical particles, energy 
and momenta conservation imply that the individual momentum at most be exchanged. 
Hence the two particle ansatz we made ought to work. 



*E.g., Tr (Z"i-iXZ"2""i-iXZ"3-"2-ixz^-"3) ^ \ni,n2,n3) 



7 



For three particles the story is radically different. This is when Integrability starts playing 
a role. Integrability means that in addition to the momentum Qi = P and energy Q2 = H, 
there exists a tower of local conserved charges Qn which commute with the momentum 
and Hamiltonian. We can introduce an arbitrary complex number u and simply encode all 
conservation laws in'^ 

00 00 

J2 [P, Qn] = 5^ [h, Qn] = . (10) 

n=0 n=0 

This relation has important consequences. For example, it allows us to guess the form of 
ip{ni,n2-,n^) in (6). The reason is that the existence of the higher conserved charges does 
imply that if we scatter three magnons with momenta {pi-,P2-,Pz\ they will scatter into some 
other momenta {^'1,^2,^3}, which must be related to the original ones by a simple reshuffling. 
In other words, the scattering is effectively pairwise. For example, for three particles we are 
thus lead to 

where . . . stand for the remaining 3 possible plane waves. The coefficient A multiplies the 
plane wave which is obtained from the plane wave with unit coefficient by swapping particles 
with momentum pi and p2- Thus 

^ = 5(p2,Pi), (12) 

where S{p,p') is the two-body S-matrix (9) derived above. The coefficient A' is the coefficient 
of the plane wave which is obtained after a sequence of two momentum exchanges, 

A' = S{p2,Pi)S{ps,Pi). (13) 

The coefficients of the three remaining plane waves are obtained in the same way. The 
generalization to > 3 particles involves A^! plane waves whose coefficients follow again the 
same pattern. Our convention for the normalization of the wave function is that the plane 
wave with no momenta exchanged has unit coefficient, i.e. 

iljim, ...,nN)= e^Pi"i+^P2"2+-+ip^n^ ^ ^(p2,pi)e^P2"i+^fi'^2+-+^P^"'^ + . . . . (14) 

Of course, this normalization depends on the choice of an ordering of the momenta; more 
on this in section 2.3. The observation that the eigenstates of the Heisenberg spin chain are 
given by such ansatz was the key insight of the seminal work of Hans Bethe [13]. States of 
the form (6), whose wave functions are given by (14) are called Bethe states. 

The energy of the multiparticle state is the sum of the energy of the individual magnons. 
Since the system is put in a finite circle of length L, the spectrum is discrete. The periodicity 
of the wave function imposes a set of quantization conditions for the N momenta, 

N 

e'P^'^l[S{p,,p,) = l (15) 



^After some n the charges are of course not independent. For the Hamihonian (7) such relation can be 
easily derived using the algebraic Bethe ansatz (ABA) formalism which we will review in subsection 2.2. 
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which are called Bethe equations. The physical meaning of this equation is the following: 
if we carry a magnon with momentum pj around the circle, the free propagation phase pjL 
plus the phase change due to the scattering with each of the other — 1 magnons must give 
a trivial phase. We shall denote by Bethe eigenstates those Bethe states whose momenta 
are quantized as in (15). Furthermore, when dealing with single trace operators, we have to 
take into account the cyclicity of the trace by imposing the zero momentum condition 

N 

Y[e'P^ = 1 



Let us stress again: the fact that the simple wave functions described above diagonalize 
the Hamiltonian (7) is absolutely remarkable and non-trivial. A generic spin chain Hamil- 
tonian will not lead to an integrable theory, the scattering will not factorize into two-body 
scattering events and the set of momenta of multiparticle states will not be conserved. Hence, 
in general, the problem will be of exponential complexity and the best we can do is diag- 
onalize small spin chains with a computer. The Bethe ansatz reduces this problem to a 
polynomial one. For example, the spectrum problem is completely solved by the simple set 
of algebraic equations (15). 

The algebraic Bethe ansatz [14], reviewed in the next subsection, provides the explanation 
for this "miracle". We will review why (10) is indeed true and we will recall that the wave 
function (14) can be constructed by acting on the "ferromagnetic" state with some beautiful 
nonlocal "creation" operators B{u). For example, the state (6) with wave function (14) is 
simply given by 

\^) = UB{u{)B{u2) . . . B{un) It ... t) (16) 

where M is some simple known normalization factor which we will write down later. The 
rapidities convenient parametrization of the momenta pj given by 

u = - cot- , e'P = '— . (17) 

2 2' u-i/2 ^ ^ 

Note that with this parametrization the Bethe equations (15) become simple polynomial 
equations 

e^^^=l where e^^. ^ ( ^^lllIlV fj ^^l^^^ , (18) 



2.2 Algebraic Bethe ansatz 

The algebraic Bethe ansatz is a formalism developed by the Leningrad school for construc- 
tively producing and solving^ integrable theories (see [10, 14, 15] for nice introductions to the 
subject and for further references). The great advantage over the coordinate Bethe ansatz is 
the constructive nature of the method and the mathematical elegance. The main drawback 

^Originally "solving" meant "computing the spectrum". 
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is that the physical picture is somehow obscured. For example, the magnon physical picture 
of (14) is somehow hidden in this formalism. 

Let us start by recalling that the Hamiltonian (7) also follows from the following object: 



^ - ^ = -8^ = *^|;('°^T.„ + >p)„ . . . (.1 + ,p).,i)__^^ . (19) 

n=l 

Before deriving this fact or explaining its relevance, let us digest the notation introduced. 
First, u is the so-called spectral parameter. It is an arbitrary complex number which we set 
to zero after taking the derivative^. The R-matrix 

Roj{u) = iuI + iF)^. (20) 

is an operator acting on a tensor product of two vector spaces: the physical spin chain 
site vector space at site j and an extra auxiliary space labeled by the index which we 
introduced. Both these spaces are isomorphic to C^, the space where the spin 1/2 lives. 
As before, P is the permutation operator that permute the spins at physical position j and 
auxiliary position 0. We can of course write the R-matrix as a simple 4x4 matrix acting on 
the vector space , |ti) , lit) All)- 

The operator inside the square brackets in (19) is the monodromy matrix 

Loiu) = Roiiu) . . . RoLiu) . (21) 

It acts on the tensor product of L + 1 spaces: the L physical spaces corresponding to the L 
lattice sites plus the auxiliary space 0. Since it is crucial to understand well the notation we 
are introducing let us be maximally pedestrian for a second. Using indices ii, . . .i^, ji, . . . ,jL 
for the physical spaces and a, b for the auxiliary space, the monodromy matrix is an object 
with indices 



The shifts by i/2 are introduced for future convenience. In other words, we can make the 
auxiliary space manifest and write 

T (v) - + V2) B{u + V2) \ 

^ C(n + z/2) V{u + z/2) ) ^^2) 

where A, B, C, V are operators which only have physical indices, i.e. they act on the physical 
spin chain Hilbert space "H = (C^)®^. 

The next object we see in (19) is the trace with respect to the auxiliary space of the 
monodromy matrix. This defines the transfer matrix 

f{u) = TroLo(M) . (23) 
^The variable u will be soon identified with u appearing in (10) and (17) in the previous section. 
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12 • • • fc ■ ■ ■ L 



Figure 2: Since i~^R{0) = P, i~^T{0) is the unit shift operator to the right, see the upper right 
corner of the figure. By definition of inverse i^T~^{0) is the unit shift operator to the left. When 
computing T'(0) the derivative will act on one of the R^s at position k, hence the sum in the last 
line. -R'(O) = I which leads to the last line. We see that T'(0) is a sum of terms which are almost 
a total shift of one unit to the right except for a small "impurity" at site k. Therefore, when 
multiplying by T~^{0), we almost get the identity operator acting on the full Hilbert space. The 
impurity simply leads to a permutation of acting on sites k and k — 1. Hence (25) leads to (7). 



Since we traced over the auxiliary space the transfer matrix is an operator acting on the 
physical Hilbert space. More explicitly we can write it in terms of the operators A and V as 

f{u)=Aiu + i/2)+V{u + i/2). (24) 

Finally, in (19) we have the derivative of the logarithm of an operator (or big matrix). This 
is understood as usual as the inverse of the matrix times its derivative. With the notation 
introduced above, (19) simply reads 



H-^, = ^,f-\0)fiO) (25) 



where the prime stands for derivative with respect to the spectral parameter u. The deriva- 
tion of (25) is illustrated in figure 2. 

An important feature of this construction is that, as explained in figure 3, transfer ma- 
trices with different spectral parameters commute, 

f{v),f{u) =0 V u,v. (26) 

The Hamiltonian is just the log derivative of the transfer matrix and hence [H,T{u)] = 0. 
Also [H,logT{u)] = 0. The latter relation is quite nice since it implies (10) neatly. Recall 
that the existence of the higher local charges Qn is what ensures Integrability of the model. 
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^02(^)^0i(m)^0i02(w - V) 




U 



t\ t\ t\ 

1 2 3 



Oi 



^0i02(w - ^^)^0i(m)^02(^) 



Figure 3: The R-matrix is a very special operator. It is designed to satisfy the Yang-Baxter 
equation depicted at the top. This equation is arguably the most important equation in quantum 
integrability. In particular it implies the LLR = RLL type relation represented at the bottom. To 
prove this relation we simply move one of the vertical lines from the left group to the right region 
using Yang-Baxter and repeat this procedure until all the vertical lines are to the right of the 
R-matrix. From this simple equation all the algebra relations (table 1) of the monodromy matrix 
elements A^B,C,T) follow trivially as explained in the main text. Finally multiplying this equation 
by and taking the trace over the tensor product of the two auxiliary spaces Oi and O2 we derive 
(26). 



Indeed, we can Taylor expand logT(u) around u = 0, 

00 

logf(^z) = 5^Q„n^ (27) 

n=0 

and generate in this way a tower of local*^ conserved charges. The first one, Qi, is the 
momentum and the second, Q2, is the energy, see figure 2. The higher terms are the higher 
charges we were looking for. Hence, indirectly, we now understand why the multi-particle 
ansatz (14) ought to work. 

Relation (26) has another interesting consequence. Since the transfer matrices commute 
with each other for different values of the spectral parameter we can diagonalize T{u) with 
a spectral parameter independent base of 2^ states such that 

f{u)\^>,) = {a{u + i/2) + V{u + V2)) 1^.) = T{uU^>,) . (28) 

^Of course we would also get conserved charges by expanding T(u) or any functional of the transfer 
matrix around any point u = u* . The advantage of expanding logr(M) around w = is that the charges 
generated in this way are local. This follows from the important relation i?(0) — P, see the derivation in 
figure 2 for an illustration of the importance of this property. 
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The states diagonalize all the higher charges at the same time since they are simple derivatives 
of the transfer matrix. We could have picked any of the higher charges as the Hamiltonian. 
The algebraic Bethe ansatz (ABA) approach tells us that all these Hamiltonians can be 
diagonalized at once with the same wave functions! We start seeing the constructive nature 
of the ABA approach. 

We already encountered the state in the coordinate Bethe ansatz approach: it is 
simply (6) with the wave function given by (14). For us it is important to construct this 
state in the ABA language. It turns out that the monodromy matrix elements B{u) can be 
used as creation operators for the magnons.^ More precisely, the claim is that 

\^) = B{u,)...BMn...'t) , (29) 

yields a state proportional to (6) once we identify Ui and Pi according to (17). In other 
words |\E') is a Bethe state. If Ui are Bethe roots satisfying the Bethe equations (18), then 
|\E') satisfies (28) and is a Bethe eigenstate. We will sometimes denote (29) by 

To understand why (28) indeed holds we need to 1) understand how to take the operators 
A and V through the creation operators B and 2) derive the action of A and V on the 
ferromagnetic vacuum. The latter is simple. From the definition (21) we easily see that 

A{u) |0) = a{u) |0) , V{u) |0) = d{u) |0) , (30) 

C{u) |0) = , B{u) |0) = a non-trivial single spin excitation, (31) 

where |0) = |t • • • t) and 
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a{u)^\^u+'-j , d{u)^(^u-'-j . (32) 

Next we need to understand the algebra of the elements of the monodromy matrix A, B, C 
and P. As explained in figure 3 the monodromy matrix L{u) satisfies the relation 

L02 (^)^oi (m)^0i02 (u-v) = i?0i02 (u - v)Loi iu)Lo^ (v) , (33) 

where Oi and O2 are two auxiliary spaces isomorphic to C^. This relation encodes the algebra 
of the monodromy matrix elements. More explicitly we can choose a basis |tt) > Iti) > lit) 7 lii) 
spanning the tensor product of these two spaces. In this basis the left-hand side of (33) simply 
reads 

fA{v) B{v) \ fA{u) B{u) \ fu-v+i \ 

A{u) B{u 
C{u) V{u) 
\ C{u) V{u)J 

while the right-hand side is obtained by multiplying the same matrices in the opposite order. 
In this way we get 16 algebra relations between the several elements which we summarize in 
table 1. 



C{v) V{v) 
Aiv) B{v\ 
\ C{v) V{v)J 



u—v i 
i u—v 

\ u-v+ij 



^This will be made more clear below, see discussion after (39). 

^'^Different overall normalizations of the R matrix (20) lead to different functions a{u) and d{u). The 
ration e{u) = a{u)/d{u) is however independent of that choice of normalization and any physical quantity 
will depend in a and d only through that ratio. With our choice of normalization for the i?-matrix we have 
B\u) = -C{u). 
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A{v)B{u 
B{v)A{u 
V{v)B{u 
B{v)V{u 
C{v)A{u 
A{v)C{u 
C{v)V{u 
V{v)C{u 

[V{v),A{u)] 
[B{u)Mv)] 

where 



f{u - v)B{u)A{v) + g{v - u)B{v)A{u) (34) 
f{u - v)A{u)B{v) + g{v - u)A{v)B{u) 

f{v - u)B{u)V{v) + g{u - v)B{v)V{u) (35) 

f{v - u)V{u)B{v) + g{u - v)V{v)B{u) 

f{v - u)A{u)C{v) + g{u - v)A{v)C{u) 

f{v - u)C{u)A{v) + g{u - v)C{v)A{u) 

f{u - v)V{u)C{v) + g{v - u)V{v)C{u) 

f{u - v)C{u)V{v) + g{v - u)C{v)V{u) 

g{u - v) [Aiv)V{u) - Aiu)V{v)] = g{u - v) [V{u)Aiv) - V{v)Aiu)] (36) 
g{u - v) [B{v)C{u) - B{u)C{v)] = g{u - v) [C{u)B{v) - C{v)B{u)] 
[Ciu),C{v)] = [A{u),Aiv)] = [V{u),V{v)] = (37) 



9[u) 



I 
u 



1 + 



u 



{3t 



Table 1: Algebra of the monodromy matrix elements which follows from the RLL = LLR 
relations described in the text and depicted in figure 26. 



To justify the identification of B and C with creation and annihilation operators, we 
compute the commutation relations of these nonlocal operators with the total spin generators 



jn=l 



a\ For that aim we consider the limit v ^ oo with u held fixed. In this limit 



we have, from the definition (21), that A{v)/a{v) ~ 1 + ^ (1 + S*^), C{v)/a{v) 
V{v)/a{v) ~ 1 + ^ (1 - 5^) and 

B{v) i 



a{v) 



S~ 



(39) 



Now, from the sixteen algebra relations in table 1 we can read for example the commutation 
relations [B{u),S^] = 2B{u) and [C{u),S^] = —2C{u) which means that B{u) is a single 
spin creation operators while C{u) is a single spin annihilation operator. We also have 
[A{u),S^] = [vlu),S^] = 0. Furthermore, 



[S+,B{u)] =A{u)-V{u) 



(40) 



Equations (39) and (40) have further important consequences which we will discuss at the 
end of the next subsection. 

Using relations (34) and (35) we learn how to carry A and V through the B operators. 
Suppose we would drop the second term in the right-hand side of (34) and (35). Then it is 
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clear that (29) would be an eigenstate of T{u — i/2) with eigenvalue 



Because of the second term in the right hand side of (34) and (35) the state |\E') is typically 
not an eigenstate of the transfer matrix. We can wonder what conditions do we need to 
impose so that the contribution of these extra terms vanishes. These conditions are nothing 
but the Bethe equations (18)!^^ 



2.3 A few comments and summary of notation 

As explained above the Bethe states (6) and (29) are proportional to each other. Let us 
write a precise relation between the two. First recall that to define the coordinate Bethe 
state normalized according to (14) we pick a particular order for the momenta. Given 
an ordering pi, . . . ,PAf for the momenta, the relation between the Bethe states normalized 
according to the coordinate Bethe ansatz (CBA) (14) and the algebraic Bethe ansatz (ABA) 
(29) isi2 

= ^{«+f 1^)^° , (42) 

where we have introduced the shorthand notation 

^ n ^("^•) ' ^^^^^"^ = n ^("^ - ^^■) (43) 

i < j 

which we will use throughout the rest of the paper. Also, from now on we should explicitly 
attach an upper index "al" or "co" where necessary to specify if we are referring to a given 
object in the algebraic or coordinate normalizations. 

Of course the normalization of the algebraic Bethe ansatz state (29) does not depend on 
the order of the momenta. After all, we see from (37) that the B{ui) operators commute 
with each other. Hence we can also use (42) to go between two coordinate Bethe state with 
different ordering of the same magnons. 

Indeed, those extra terms would lead to terms of the form X]fcLi(Q^fe+'^fc) ^i"^) Y\j^k '^(^i) It • • ■ t) where 
Ofc is the contribution coming from A{u) and 5k is the contribution coming from 'D{u) in (28). Suppose we 
want to find ak- We start with A{u)B{uk) Wj^k ^i'^i) It • ■ • t) and need to carry A{u) to the right till it 
hits the vacuum. Note that we ordered the creation operators so that B{uk) comes first. We can always do 
it since [B{u),B{v)\ = 0. Now it is clear that to end up with terms contributing to we must first use the 
second term in (34) to commute A{u) through B{uk) so that B{uk) becomes B{u) {A{u) becomes A{uk)). But 
then we must always use the first term in (34) to commute A{uk) through all other S's since we are already 
happy with the arguments of the B operators. Hence we conclude that = a{uk)g{u — Uk) Ilj^/c fi^j ~''^k)- 
Similarly 5k = d{uk)g{uk — u) Ylf^k fi^k — Uj). The condition ak + 5k = yields the Bethe equations (18) 
once we recall (38) and (32). 

^^In the local basis of states (6), the wave functions are rational function of the rapidities {ui}. The 
conversion factor in (42) can be simply derived by demanding that the two wave functions have the same 
zeros, poles and large Uk behavior. 
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Let us re-write the Konishi state (4) in the coordinate and algebraic normahzations. The 
Konishi state has L = A, N = 2 and momenta pi = —p2 = 27r/3 or mi = —U2 = y/S/G. This 
leads to 



K = 4e-^^/^ 



Tt{ZZXX) - Tt{ZXZX) 



K = — 
27 



Tt{ZZXX) - Tt{ZXZX) 



in the CBA normalization, 
in the ABA normalization. 



Finally let us end this section with a discussion on the completeness of the basis of states 
discussed so far. Using (40) and the same kind of reasonings of footnote 11 we can easily 
prove that 

S+Biu,)...BMn...t) = (44) 

if the roots uj obey Bethe equations. In other words, kills Bethe eigenstates. This means 
that Bethe eigenstates are highest weight states. Indeed, if we count the solutions to the 
Bethe equations (18) we find that there are precisely [14] 



N \N-1 



solutions with spin flips. This is exactly the number of highest weights for a chain of 
length L and = L/2 — N. All other states are found by acting with on these states. 
Indeed 

L/2 I- 



E 

N=0 

Hence a complete basis of states is given by 



2^ . (45) 



{S-y \{u,}r = (5-)" S(«0 . . . BM It ... t) (46) 

where Uj obey Bethe equations and n = 0, . . . ,L — 2N. All these states have the same 
energy (or any other charge) as the highest weight state B{ui) . . . B{um) It ... t) since the 
Hamiltonian (or the transfer matrix) commutes with for i = z,+,—. 

Note also that these states can be written using only the creation operators B{u) since 
for large rapidities this operator becomes the lowering operator, see (39). More precisely, 
this state is simply given by 

{S-y \{u,}r' = hm I n ^ ) B{A^) . . . B{K)B{u,) . . . B^ |t . . . t) . (47) 



In short, we see that if we consider the Bethe equations for uj and then add a few Bethe 
roots at inflnity we describe the full Hilbert space. 

The factor (11^=1 ^l"^/"^) simply ensures a good limit when the roots go to inflnity. The 
coordinate Bethe ansatz states have a better limit when we send some of the roots u^: to 
inflnity 

hm \{u,}r = S-\{u,},^,r- (48) 
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Particularly important states are the so-called vacuum descendants 

|{oo^}r = (5-)^|t...t) (49) 

which correspond to a state with roots at infinity only. They correspond to a ferromagnetic 
vacuum rotated away from the z axis. In the = 4 language it corresponds to an operator 
with L — N scalar fields Z and scalar fields X without any anomalous dimension. These 
are BPS states whose anomalous dimension is zero to all orders in perturbation theory as a 
consequence of supersymmetry. In the string theory dual language these states correspond 
to BMN point like strings which rotate around one of the equators of at the speed of light. 
Different values of correspond to different equators of which are of course related by a 
global rotation. 

3 Tailoring tools 

In this section we will present the necessary tools for computing the three-point function 
in the fashion explained in the introduction. That is, we will express a closed chain Bethe 
eigenstate as an entangled state in the tensor product of the two open subchains Hilbert 
spaces (subsection 3.1). Then we will explain the flipping procedure, mapping a state in one 
of the two subchains from ket to bra (subsection 3.2). Finally, we will compute the overlap 
of Bethe states (subsection 3.3). We denote these three procedures by cutting, flipping and 
sewing. In section 4 these tools will be used to put the pieces together into the structure 
constant as illustrated in flgure 1. Two references that give a nice, pedagogical introduction 
to the techniques used in the cutting and sewing sections below (from the algebraic basis 
point of view) are [10, 15]. 

3.1 Cutting... 

Consider a spin chain of length L and a generic state (6) of that spin chain. We denote the 
first I spins starting from the left of that chain by lefl subchain and the last r = L — I spins 
by right subchain. We can represent that state in the big chain (6) as an entangled state in 
the tensor product of the left and right subchains as 

rain{N,l} 

1^) = ^ ^ ^ ^(rii, ... ,raAr) 1^1,^2, ... ,rifc) (g) l^fc+i - Z, ... ,raAr -/) . 

A;=0 l<ni<---<ni.<n<r!,/;_|.i<---<n]v 

(50) 

The sum over k is the sum over how many magnons are in the left chain. The sum on 
the right hand side of (50) has (j^) terms. For Bethe states there is a huge simplification. 
Namely, when we represent a Bethe state as an entangled state in the two subchains, each 
of the subchain states still has the same Bethe state form. 

The reason is that - locally - these are eigenstates of a local Hamiltonian and therefore 
take the specific local form presented in the previous section. As a result, when we represent 
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a Bethe state as an entangled state in the two subchains, each of the subchain states still 
has the same local Bethe state form. 

In other words, a magnon that is locally propagating along the chain does not know 
that at some far away point the chain was broken. To write the corresponding piece of the 
wave function, all we need to know is whether that magnon propagates on the left or right 
subchains. A Bethe state therefore breaks into two as 

\{u.i}) ='^H{a,a)\a)i'S)\a)r , (51) 

where the sum is over all 2^ possible way of splitting the rapidities into two groups a and 
a such that aU a = {ui}. For example, if = 2, the possible partitions {a, a) would be 
({}, U2}) , {{ui}, {U2}) , {{U2}, {ui}) , U2}, {}). As opposed to the (^) terms in (50), 
we only have 2-^ terms in (51). This is a very convenient simplification whenever N . 

The function H{a,a) takes different forms in the coordinate and algebraic bases. Below 
we will compute that function for each of these normalizations. Before moving on, let us 
introduce some further notation in the same spirit of (43). We shall use 

^ n ^(«^) ' F-^^W F{u. - V,) , = n ^("^ - • (52) 

u-j&oi Ui a Ui,Uj a a 

Vj £ a i < j 



Cutting a coordinate Bethe state 

The normalization of a Bethe state in the coordinate base depends on the choice 

of ordering ui,U2, ■ ■ ■ ,un of the magnons, see (14). For the states in the two subchains 
and we choose the ordering induced from the ordering of With these 

conventions we will now derive H{a, a). 

We first consider the case where all the magnons are in the left subchain: a = {ui}, 
a = 0. In this case, by construction, the wave function of coincides with the wave 
function of and therefore 

if^°(K},0) = l. 

Now suppose we shift some set of magnons from a to a. There will be two factors contributing 
to H. One is the phase shift acquired by the a magnons when translated trough the first I 
sites^'^ 

where ai (a^) and di (dr) are defined as in (32) but using the lengths of the left (right) 
subchain instead of L. The second contribution to H{a, a) is the scattering phase between 

-"^^In other words, this factor arises because we shifted the label of the positions of the magnons of the right 
chain by I in order to start counting them from 1 as usual, see (50). 
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all pairs of magnons Uj G a and G a such that i < j 

n jT^y (53) 

J- J- f(^Ui-Uj) 

I < J 
Uj G a, Ui G a 

Multiplying the two factors and using the shorthand notation introduced in (52), we can 
write H'^'^{a, a) as 

„cS £aa faa faa 

H-^{a,a) = % ^ i</< . (54) 
Finally, note that cutting a descendant of a coordinate state is trivially obtained from (48) 

\{ui, oo"})™ = J2 ^"(«' J2 Q I" ^ ® I" U {oo}"-'")™ . (55) 



Cutting an algebraic Bethe state 



We can now use the translation between the coordinate and the algebraic normalizations 
(42) to conclude that 

H''\a,a) = r^'dyf . (56) 

We can also straightforwardly derive (56) by writing the monodromy matrix (21) as a product 
of the monodromy matrices of the left and right subchains, L = (_Roi • • • Roi) {Roi+i ■ ■ ■ Rol) = 
LiLr . More explicitly, the monodromy matrix elements are split as 

fAiu) Biu)\ _ fMu) Biiu)\ fAriu) Briu)\ 

[Ciu) V{U)) \Cl{u) Vi{u)) \Cr{u) Vriu)] ' ' ' 

In particular, the Bethe state (29) can be written as 

N 

^(u)|0) = n {Ai{uj)Br{u,) + Bi{u,)Vr{u,)) |0) . (58) 

We could open the parentheses and get a sum of 2^ terms parametrized by which roots end 
up in the left spin chain creation operator Bi. That set of roots is denoted by a. The set 
of roots ending up in the right spin chain creation operator Br is denoted by a. Of course 
{uj} = a U a. This is precisely the sum over partitions in (51). To get rid of the A and 
T> operator in (58) we use the commutation relations in table 1 to commute them through 
the i3's plus the actions of these operators on the vacuum of the corresponding chains as 
given in (31). In this way one arrives at (56). This technique is known as the generalized 
two-component model [16, 41]. 
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3.2 ... Flipping 



After splitting the states into two (which we did in the previous subsection), we need to 
represent the Wick contraction of the elementary fields as a spin chain operation (see figure 
1). The Wick contraction operation takes two ket states on two subchains and produce a 
number. It can be achieved by two steps. Wc first flip one of the two kcts into bra and 
then, in the next section, contract the ket and bra states. The result does not depend on 
which of the two kets we choose to flip. In what follows we will always flip the state in the 
right subchain before contracting it with the state on the left subchain of the next operator 
(as illustrated in figure lb). As explained in the introduction, the flipping procedure is 
not the usual conjugation. The usual conjugation acts on the corresponding operator by 
conjugation and therefore fiips the order of fields and their charges. On the contrary the 
hipping procedure J-' does not change the operator. As a result, the contraction of a ket state 
with a flipped bra state is the same as the Wick contraction of the corresponding operators. 

In the introduction we added an upper arrow to distinguish ;(\I' | = To |^); from ;(^'| = 
(1^);)^ For example, for real we have 

t : e''''\XZXZZ) ^ {XZXZZ\e-''^ , (59) 

T : e''*'\XZXZZ) ^ {ZZXZX\e^'^ . (60) 

Note that in our notations {Z\Z) = {Z\Z) = 1, {Z\Z) = and the spins in the kets and bras 
are ordered from left to right, i.e. {ZX\ZX) — 1 while {XZ\ZX) — 0. Summarizing, the 
difference between the two maps is most clearly illustrated by their action on the local basis 

t : ^(ni, . . . ,n7v) |ni, . . . ,nAr) H- ^''^('^i, • • • , ^iv) (^i, • • • , ?^jv| , 
T : i){ni, . . . ,nN)\ni, . . . ^un) ^ i){ni, . . . ^un) {L — un + 1, ■ ■ ■ ,L — ni ^-1\C , 

where C stands for charge conjugation, which exchanges Z -k^ Z and X -h^ X. 

We can now flip any Bethe state. Since the Hamiltonian is invariant under flipping the 
orientation of the chain and the charges, the operation T maps a ket Bethe state into another 
bra Bethe state. The overall factor depends on the normalization. To determine that overall 
factor in the coordinate normalization, consider flrst a two-magnon state 

U2}Y° = [e'^'i"^+'^'"' + ^(P2,pi)e'^^"^+'^^"^] |ni, n^) . 

By changing variables to mi = L — n2 + 1 and rn-i = L — ni + 1 it is easy to see that 

This construction is trivially generalized to any number of magnons. We find 

jTo \{ui\Y° = e}"^ . > , , "°{{<}\C. 
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Our goal was to flip the second ket in (51) which we can now do. We have 

„Q:— I £aa £aa £aa 

•F^econdchainOlKir = E ^2 WW I ^ ^° (c^l ^ • (61) 

aUae{u} 3 ^ /< 

where a* is the set a with its elements complex conjugated. Similar expressions can be 
written in the algebraic conventions using the conversion factors derived before. ^'^ 



3.3 ... and Sewing 

The last building block which we need to understand is the overlap of Bethe wave functions, 
i.e. scalar products. The procedure of overlapping the wave functions is what we denote as 
sewing. The importance of spin chain scalar products in = 4 SYM was pointed out in [4]. 

The quantity we are interested in is^^ 

N N 

Sn{{v},{u}) ^ {0\l[C{v,)l[l3{u,)\0) . (64) 

i=i i=i 

The explicit form of this quantity is actually known [17]. It is written for completeness in 
(114) in appendix A. If {u}, {v} take particular values this scalar product simplifies dra- 
matically (luckily!). For example, when we consider the norm of a Bethe eigenstate we 
have {u} = {v} and the Bethe roots obey the Bethe equations. In this case we have the 
remarkably simple result [18] 

5^({n}, {u}) = rf^«> a^">/i"»"Vi"^^"^ det 9,0, (65) 

where 0, was introduced in (18), dj = and we are using the shorthand notation introduced 
in (43). We shall also use detdjcpk = det0'r i. In the coordinate basis this formula reads 

1 /■{«}{"} 

U-({u}) . = det0;„, . (66) 

As explained in section 2, descendants of Bethe states are obtained by sending some of the 
roots to infinity (48). We use to denote the total number of roots and M to denote the 



^^Using (42) we see that the algebraic Bethe states transforms as 

J-o|M)-i = (-l)^-i(K}|C. (62) 

Furthermore 

J-second chain O |M)^l = T " rf>f 1 0) , $5. (0 |C"* (7 . (63) 

aUa(z{u} 

Recall that in our normalization C{u*) = — [B{u)]^ . 

^Ht is trivially related to the scalar product through S'^{{v*}, {u}) = (—1)^ ^H{v}\{u})^^ which in turn 
follows from C{u*) = - [-B(u)]^ 
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number of finite roots. The corresponding norm of a Bethe eigenstate descendant is related 
to (66) by the square of the Clebsch-Gordan coefficients 

co/r„ ^^N-M^\r„ ^^JV-A/i\co ("^ ~ ^^)' ~ co/ r„ i i r„ i \co /«^7^l 

{{u;oo }\{u;oo |) = (^i ^ m - N)\ u^jlW/ • ^^') 

Another notable simplification occurs when we send one set of Bethe roots, e.g. the vj, 
to infinity. In this case we find^*^ 

-({oo}|{^})-= ^I^./mW E (-l)'"'e"r^ (68) 

3 J < a\Ja={u} 

where the sum runs over all possible partitions of the N elements of the set {u} into subsets 
a,d and \a\ denotes the number of elements in a. Similarly, the inner product of a Bethe 
state descendant with a vacuum descendant is related to (68) by a product of Clebsch-Gordan 
coefficients 

-({oo}|{^, oo^-^})- = ^°({oo}|{n})- . (69) 

In this formula the number of roots u is M. 

There is one more case where the general formula (64) simplifies. This is the case when 
one of the operators is a Bethe eigenstate while the other one is left generic. This is presented 
in (119) in appendix A. 

Note that sewing Bethe states involves at most sums over partitions of the magnons. 
The number of terms therefore only grows as a power of the number of magnons and is 
independent of the chains' lengths. 



A new recursion relation for scalar products 

For completeness let us illustrate how we could compute easily (64) for arbitrary {u}, {v}. 
One possible strategy for computing this quantity is to derive a recursion relation yielding 
Sn in terms of Sn~i- Such recursion relations exist in the literature, see e.g. [17], and are 
presented in appendix A for completeness. 

Here we derive a new recursion relation which we found quite efficient computationally. 
The idea is to get rid of particle creation operator B{ui) by writing 



Sn{{v},M) = {0\ 



N 



l[C{v,),13{u,) 



N 



l[B{u,)\0) 



(70) 



^^Note that in this case the algebraic Bethe states are not well normalized. Indeed, as we saw before, 
Ciuj) ~ ia{uj)/uj S'^ — >■ c» as Uj — oo. The coordinate wave functions are better normalized and do not 
vanish or diverge in this limit (48). 
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which holds since the B operator kills the bra vacuum. Next we use the algebra relations 
(36) to compute this commutator. We have 



■ N 



(n-l \ / N 

l[C{Vj)\g{Ui-Vn){A{Vn)V{m)-A{u,)V{Vn))i H 
j=l / \j=n+l 

(71) 

The term between the two products is just [C{vn),B{ui)]. Now we pick the operators A and 
V in this expression and carry them to the left through the Cs using the algebra relations. 
Recall that A and V have a trivial action on the left vacuum. This procedure does not change 
the number of Cs which is now equal to — 1. So, what are the possible final arguments of 
the remaining iV — 1 Cs? When commuting an A through a C we have 

C{v)A{u) = f{v - u)A{u)C{v) + g{u - v)A{v)C{u) , (72) 

which means that the A can keep the same argument or swap arguments with the C The 
same is true for a V operator passing through a C It is clear that we will get terms where 
the arguments of the Cs are — 1 out of the N original f 's. There is however one more 
possibility which comes from the algebra terms which swap the argument of the operators: 
we can also end up with — 2 of the original u's together with the root ui. That is 

Sn {{Vl, Vn}, {Ul, Un}) = ^K Sn-i{{Vi, ...,Vn,..., Vn}, {Ul, U2,..., Un}) 

n 

- ^ Cn^m Sn-i{{Ui,Vi, . . . Vn, ■ ■ ■ , Vm, ■ ■ ■ Vn}, {^1, ^2, • • • , Un}), (73) 
n<m 

where the hat means that the corresponding root is omitted. Naturally 5*0 = 1. We simply 
need to find 6„ and c„,m! Let us derive 6„. We start by re-ordering the Cs in a smart way 
and write (71) as^^, 

l[C{v,)\ C{vn),B{u,) 

We are interested in 6„. In other words we want to consider the contribution to the scalar 
product where C(w„) and C{ui) are not present. This means that we must get rid of C(f„) 
in (74) and therefore the only term which will contribute is 

nC(t;,) J [C{vn),B{ui)] = \ \{C{vM g{ui-vn) {A{vn)V{u,) - A{Ui)V{Vn)) . (75) 

Now the ^'s must travel to the left using (72). Furthermore, when using these relations we 
must always pick the first term in the right-hand side. This is the term where the A and C 
do not swap their arguments. Otherwise we would end up with a C with an argument f„ or 
Ul and therefore this would not contribute to 6„. For example, 

N / ^ \ 

AM=a{vn)l[f{v,-Vn){0\ n'^(^^)) +•■■ (^6) 
^''Note that we can order them in any way we want since they commute, [C{u),C{v)] = 0, see table 1. 



(74) 
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where . . . stand for terms where one of the Cs ended up with the argument f„. We used 
the fact that the action of A{vn) on the bra vacuum simply yields a(v„). The other A and 
V operators in (75) are treated in the same way. It should be clear that at the end of the 
day we get 

bn = g{ui - Vn)a{Vn)d{Ui) J]^ f {Ui - Vj)f{Vj - Vn) + («! O f„) , (77) 

With a similar reasoning we could derive c„,m. We would find 

N 

Cn,m = 9{Ul - Vn) g{Ui - Vm) a{v^)d{Vn) f {Vn " t^m) JJ f{Vn - Vj) f{Vj - Vm) + {u ^ Ul) .(78) 

Since we know how to convert between the coordinate and algebraic Bethe states we can easily 
write a recursion relation for scalar products in the coordinate Bethe ansatz normalization. 
This is presented in appendix A. 

The recursion relation (73) with (77) and (78) provides a complete solution to any scalar 
product in a straightforward way. For example, it is straightforward to implement this 
recursion in Mathematica. As mentioned above, there exists another recursion for the general 
scalar product (64) in the literature, which we present for completeness in (118) in appendix 
A. 

It is interesting to compare the two recursions in the following table: 



New recursion (73) 


Usual recursion (118) 


Uses scalar products with less particles as 
fundamental building blocks. 


Uses generalized objects as building blocks. 
I.e., scalar products in different theories with 
different a and d functions appear at every 
step of the recursion and less particles. 


Derived by reducing the number of particles 
in a very explicit way. Notion of particle is 
fundamental. 


Derived from the analytic properties of the 
result. The scalar product is a rational func- 
tion which can be therefore reconstructed 
from its poles and zeros. The notion of re- 
ducing the number of particles is secondary. 


Contains much less terms than a brute force 
computation of the scalar product but still 
more terms than the recursion (118). 


Contains much less terms than a brute force 
computation of the scalar product and less 
terms than the recursion (73). 



Given this comparison it is fun to draw an analogy with two very important recursion 
relations in = 4 SYM which are used in studying scattering amplitudes: the CSW recur- 
sion relations [19] and the BCFW expansion [20]. A very similar table would be suitable for 
comparing these two recursions. The CSW and BCFW recursions would be the analogue of 
(73) and (118) respectively. 
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4 Three-point functions 



Having introduced the necessary tools in the previous section, we are now ready to put 
them together to compute the objects of interest of this paper: structure constants C123 of 
three single trace operators in A/" = 4 SYM at tree level. An alternative way or writing the 
structure constants in a scheme independent way is as 



C 



123 



(79) 



where the indices are identified modulo 3. The right hand side is indeed a Xj independent 
quantity which coincides with C123 once we use (1) and (2). We want to study these structure 
constants at tree level, i.e. we are after the first term in the expansion (3). We will focus on 
a particular subset of operators as we will now describe. 



4.1 Setup 

The operators we will consider are linear combinations of single trace operators made out 
of two complex scalars. As reviewed in the previous sections, such operators with definite 
one-loop anomalous dimensions can be represented by Bethe eigenstates on a spin ^ chain. 
This fact will allows us to use the technology introduced in the previous section for cutting 
and sewing the operators. 

The general setup that we will use in the rest of the paper is presented in figure 4a. We 
consider three operators O2 and Os, with corresponding lengths Li, L2 and L3. The 
number of propagators between operators i and j is fixed to be 

lij = ]^{U + L,-Lk). (80) 

where (i, j, k) is a permutation of (1, 2, 3). We will restrict ourselves to the non-extremal case, 
that is where all /jj's are strictly positive. This is not only the generic case but it also has one 
important advantage: for non-extremal correlators the contribution of the operator mixing 
with double trace operators is suppressed in 1/A^c and need not be considered. Furthermore, 
there is considerable evidence in the literature that the extremal correlators can be considered 
as an analytic continuation of the non-extremal ones when some lij — )■ [21]. 

Each operator Oi is a single trace operator made out of two complex scalars and is mapped 
to a spin chain state. One of the scalars we interpret as being the vacuum (or spin up) while 
the other scalar we interpret as excitations on that vacuum (or spin down). The total length 
of the operator is Li and the number of excitations is denoted by Ni. The different scalars 
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Z X 

02 



a 



{v} 
(b) 



Figure 4: (a) Three-point function of SU{2) operators at tree level. All contractions are 
such that R-charge is preserved. This is the simplest non-trivial configuration which is not 
extremal. Note that the number of excitations on each chain is subject to the condition 
Ni = N2 + N3. Also, if we denote by lij the number of propagators between operators i and 
j, we have I12 = Li — N3, /13 = N3 and I23 = L3 — N3. (b) We show the partitions of the 
excitations in Oi. Note that we only need to use the operation of cutting for this operator, 
as a part of its excitations (a) are contracted with those of O2 and another part (a) with 
those of O3. 



for the different operators are summarized in the following table: 

vacuum excitations 





z 


X 


02 


z 


X 


03 


z 


X 



(81) 



Note that this is the only setup that is fully contained in the SU{2) sector and which involves 
non-extremal correlators at the same time. When Wick contracting the single trace operators 
with each other we can only connect a scalar with its conjugate. Therefore we see that all 
the vacuum constituents of are connected O2 while all excitations of O3 are contracted 
with Oi, i.e. 

l23 = L3-N3 , /i3 = Ar3. (82) 
The richest contraction is between operators Oi and O2 at length 

h2 = L,-N3. (83) 

-•^^Note that there is no physical difference between Oi and 02- Indeed, our results will be invariant under 
the exchange of the two. The choice of vacuum (81) however is not symmetric. Still, we found this choice 
more convenient to work with than a more symmetric one. 
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In this case we can overlap both vacuum and excitations. Finally, note that the total number 
of excitations in the three operators are constrained by 

Ni = N2 + N3. (84) 

These statements are illustrated in figure 4a. 



4.2 Brute force computation 

Representing each operator by a spin chain state of the form (6), we can always compute 
the three-point function of figure 4a by brute force as 

Cl23 = '^il..mN2M-N3+l,Li-N3+2,...M'^L2+l-nN^,...,L2+l-ni'^l,2,...,N3 (§5) 

where 

Li L2 L3 

a 



l<ni<-<nN^<Li-N3 



and tp^^^ is the wave function of Oj given in (14). The factor a contains some symmetry 
factors and the normalization of the states N'j which is given by 

A/; = E ■ (86) 

l<7ii<---<njv^ <Lj 

This way of computing the structure constant is equivalent to the direct Feynman diagram 
computation and involves the precise form of the iV-particle wave function (14), which has 
A^! plane wave terms, and a sum of O ((Li — N^^)'^^^ terms. For large number of excitations 
and long chains, (85) is very inefficient. This is when the techniques we introduced in the 
previous section will prove most powerful. Namely, we will see that using the operations of 
cutting, flipping and sewing, one can express the structure constant purely in terms of the 
rapidities and lengths of each spin chain. 

Note that the absolute value of the structure constants is independent of the normalization 
of the three operators since we divide the overlap by the normalization of the wave functions. 
On the other hand changing the phase of these wave functions does not change their norm 
but it will modify the overlap by a phase. Therefore the phase of the structure constants 
is sensitive to the normalization of the operators. The phase is important however for 
bootstrapping higher-point functions and therefore we will keep track of it.^^ To fix this 
phase we chose to work with the coordinate convention for the wave functions. 



^^An analogy is the metric on curved spacetime which depends on the choice of coordinate system. 

^"Recah that the coordinate normahzation is sensitive to the order of Bethe roots. For different orders we 
get therefore different structure constants. Of course, the absolute value of C123 is always the same but the 
phase does change. 
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4.3 Tailoring three-point functions 



To compute the three-point function using the tools introduced in section 3, we first need 
to use the cutting and flipping operations to cut the three operators and prepare them to 
be sewed. In principle each operator should be decomposed and flipped as in (61). For the 
operator Oi the lengths of the left and right subchains in this formula are k-i^i and 
respectively (80), (where the indices are identified modulo 3). Recall that in this formula a 
is the subset of magnons on the left subchain and a is the subset of magnons on the right 
subchain. Now, as depicted in figure 4, only very particular subsets contribute when we cut 
the operators and 02- Namely we need as = for and a2 = for 02- The non-trivial 
cutting is that of Oi. The length of the left and right spin chains in this case are Li — N3 
and respectively. 

The second step is to compute the scalar products between the different subchains, this 
is what we called the sewing procedure above. Note that the contractions between operators 
O2 and O3 are trivial, as we are simply contracting vacuum fields. On the other hand, the 
contractions between Oi and O2 and between Oi and O3 are nontrivial and we need to 
use the inner products of Bethe states to compute them. Finally, we normalize our result 
dividing it by the norm jVj of each Oj. We use the notation shown in figure 4b and denote 
the rapidities of Ci, O2, O3 by u, v, w respectively. 

At the end of the day, the three-point function is given by 

''^ = f<Kwl> E <«/"/>"/r({«-}i">(<s-i{»}> (87) 

•I < J < QUa={n} 

where ^ 

di[u) \u-lJ u 

In these formulae all quantities are computed in the coordinate normalization as mentioned 
above. We chose to omit the explicit "co" superscript to make the formula less cluttered. 
So, 

(«1 W) = W)^3 ' = L,-^:(K}l«)i°-^3 • (89) 

and 

M = i>|«)i° , M2 = Z(v\v)z , ^f3 = Z{Mw)Z (90) 

Our expression (87) is completely given in terms of inner products between Bethe states 
in the coordinate normalization. These inner products can be found using the recursion 
relation derived in the previous section or using the general formula (115). Thus we solved 
the problem of computing c^23 ^'^^ generic states in our setup. 

^^As explained above, different normalizations of the wave functions yield the same absolute value for C123 
but different phase for the structure constant. In other words, we get the same absolute value for C123 for 
the coordinate wave function with any order of Bethe roots or even for states in the algebraic normalization. 
But the phase will differ for all these cases. 
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Furthermore, most inner products in (87) are not the most generic ones and therefore 
the formula can be simphfied considerably. For example, for the normalization factors Mj 
we can simply use (66) which is of course much simpler than the generic inner product! The 
inner product {q.*\{w}) is also quite simple: since the states \{w}) and \q*) are proportional 
the vacuum descandant with all spins down (see figure 4b), the product {q.*\{w}) factorizes 
into^^ 

{a*\{oor-^){{oor-^\{w}) 
"^"^^ = ({oo}-3|{oo}-3) • (91) 

For the inner products in this expression we can now use (68) and (69) which are again much 
simpler than the general case (115). In sum, the only complicated inner product is ({w }|a). 
For this case we should indeed use (115) generically. There are some cases of course when 
even this inner product simplifies. For example, if we send the roots v to infinity or if we 
send the roots u (and thus a) to infinity. In other words, when either Oi or O2 is a BPS 
operator. 



4.4 General cases 

In this section we summarize the final results for generic operators Oi. To present the results 
in a concise way we should introduce some useful notation. As before we use Li to denote 
the length of operator Oj, Ni to denote the number of excitations of this operator and Mj to 
denote the number of finite Bethe roots. The sets of finite roots of the three operators are 
denoted by {«}, {f }, {w}. 



We further introduce 



^Jg^-^g^-^\ f> '''''' 



(92) 



^(/|M) ^ J2 i-lPrVe? (93) 

aUa={u} 

where ipj is introduced in (18), ai,di,f are defined in (88) and a and a are partitions of 
{u}. As usual, we are using the shorthand notation (43) and (52). The set of Bethe roots 
appearing in is for a = 1,2,3 respectively; furthermore recall that these 

sets comprise only the Ma finite roots of each operator. Finally we denote 

^(/|{'u},{t;}) = Yl 9<^^''^9>^^''\-lf'^''''e^e'lh^''h''^h^^h''''dett^''detr^ . 

aU a = {v} 

where (— 1)'^" is defined as a sign of permutation of the ordered set {v} which gives a U a^^. 
Finally 

h(u) = l-iu , t(u) = — r (94) 

U(U + l) 



Again we omit the subscript and the superscript "co" which is common to all ket's and bra's in this 
expression. 

^■^In Mathematica we define (— 1)"^° as sign[a_List, ab_List, v_List] := Signature[Join[a, ab]]Signature[v] 
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yOOOO 

•-^123 



'-^123 



'-^123 



1 



'/l2 

N2 



N2 



M2 
M2 



^{h2\{v})^{h3\{w}) 



'-^123 



'-^123 



"-^123 



H2 — M2\ ^{lulM) 

N2-M2J ^i^2^3 
'■12 ~ 



No 



■>li^2^3 



h2-M^\ ^(/i3|W)^(/l3|W) 

N2 



?1^»2^3 



= ^ 4T"^ih3\aWih2\a, {v}) 

aU a = {u} 

aU a = {u} 
l"| = /l3 

Table 2: Tree level structure constants 0^23 for three SU{2) operators in the setup of figure 
4. We use o to indicate a BPS state and • to label a non-BPS state. So, for example, C^** 
corresponds to the structure constants when Oi is a protected BPS operator while O2 and 
O3 are generic Bethe eigenstates. The most general case is of course C"*- The last two cases 
were computed for highest weights only, i.e. for Ma = Na- It is straightforward to generalize 
them to the case comprising descendants but the formulas become more involved. The first 
two cases are discussed in greater detail in the next subsection. 



In addition to our usual shorthand notation we are also using 

dett^"= det t{ui-Vj) (95) 

e a 

In terms of these useful functions the final results for 0^23 are presented in table 2. Note that 
C123 and (7*23 as well as between C^'s C*23 are not different cases. That is, Oi and O2 
are exchanged if we choose to view the same operators as Z's {Z^s) excitation of the X's 
(X's) instead of X's (X's) excitations of the Z's (Z's). We notice that C°23 and C^** 
only nonzero when M2 < I12 < -^2 — ^2 due to cancelations in ^2(^12)- 

We now move to the discussion of the two simplest cases which are the cases when all 
operators are BPS or when only O2 is not BPS. 



4.5 Simplest examples 

The simplest three-point functions involves three protected (BPS) operators. This case is 
protected by supersymmetry and was studied in [22]. The tree-level result holds at any 
coupling. In our language this correlation function involves three operators where all exci- 
tations have zero momentum, i.e. all Bethe roots are sent to infinity. In other words, all 
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three operators are vacuum descendants in the SU{2) sense. In this case (87) simphfies to a 
simple combinatorial factor^^ 

Li-N- 



\J L1L2L3 



cfl = ^^123 = \ ^ . (96) 

-^1 \ / -^2 \ / -^3 



The next to the simplest case is when one of the operators is not protected. In this case 
the three-point function is no longer fixed by (super) symmetry and it has a very interesting 
and rich structure. Consider for example the case when Oi and are BPS states (i.e. 
vacuum descendants) while O2 is a generic Bethe eigenstate. In this case, the structure 
constant c^23 be a function of the Bethe roots characterizing the operator 02- If all 
these roots are finite the operator corresponds to an highest weight spin chain state. When 
N2 — M2 of these roots are sent to infinity we obtain SU (2) descendants. I.e. we generate 
the full SU{2) multiplets by acting on the highest weights with (^g-^^2-M2 ^ rpj^^ number of 
finite roots is M2 and we use {v} to denote the set of finite roots. The structure constant 
0^23 = C^'s for the case when only O2 is non-BPS is then given by 



It T T (Li-Nz-M2\ 
^0.0 _ V^1^2^3i ^2-M2 ) 
'-"123 



■N3 



(97) 



There are several interesting limits which we can take in (97). We can consider many 
magnons, few magnons, the near BMN limit, the classical limit etc. 

For example, the simplest possible nontrivial case is when O2 has only two particles with 
opposite momenta. In this case we have simply M2 = N2 = 2 and vi = ~V2 = v = \ cot |. 
Furthermore, from the Bethe equations (15) we can see that in this case p = 27m/{L2 — 1), 
with n G Z. Then 



^0.0 _ -W2 / L,L2L, sin(^) sin(|(/,2-l)) 

p{m)o -Ml) ' 

where lu = {Li + L2 - L^)/2, N3 = {Li + L3 - L2)/2 and Ni = N3 + 2.^^ 

The BMN limit is discussed in appendix B. We briefly comment on the classical limit 
in the next subsection. A more thorough analysis, comprising in particular a thorough 
discussion of the classical limit of all our results, will be reported elsewhere [23]. 

^^The normalization of the coordinate Bethe ansatz state when all roots are sent to infinity is given by 
™({oo}^|{oo}~)™ = see (69). 

^^Note that (98) is invariant under I12 —> I32 — {L2 + Li — Li)/2 which is a particle hole symmetry. 
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4.6 Scaling limit 



The scaling (or thermodynamical) limit is also known as the Sutherland limit [24] and corre- 
sponds to low lying excitations around the ferromagnetic vacuum. In the AdS/CFT context 
this limit was re-discovered in [25] and proves to be a very useful limit for comparison with 
string theory computations through the so-called Frolov-Tseytlin limit [26] which we describe 
below. In the scaling limit the length L and the number of Bethe roots M is very large. The 
Bethe roots also scale to infinity as 

M ~ L ~ ^ oo (99) 

and are of the same order. In this limit the roots distribute themselves in umbrella shaped 
contours described by some density p(f ), see figure 5. These contours are also called cuts. 
A systematic description of the spectrum in the classical limit is achieved through the finite 
gap method of KMMZ [27]. In this method one introduces a resolvent 

G{u) = y ^— = [ dv^ (100) 

where the integral is over all the contours Ck where the roots lie in the continuum limit. The 
(derivative of) this function describes a hyperelliptic curve. Our result can be written very 
neatly in this language. We find that the structure constants reduce to nice integrals on this 
Riemann surface, [23]. For example, for the structure constants C°°*, C°*° and C*°° we find 
that the building blocks 

:^^~exp dt i j ^^q\og{l - e^^') - j ci«p log (2 sinh(7rtp)) | (101) 
(up to a phase factor) where^^ 

qiu) = -- Giu) (102) 
u 

is a sort of trimmed quasi-momenta, (see figure 5). In this expression all the roots are 
taken to be finite. Expression (101) holds for any scaling configuration of roots {v} with 
any number of cuts.^*^ The scaling limit of the more general cases in table 2 can be also 
analyzed. A detailed analysis of the classical limit and other interesting regimes will be 
published elsewhere. 

For concreteness let us consider the case where O2 is non BPS while Oi and are 
protected. The simplest example of the configuration with the scaling eq.(99) is the con- 
figuration with one cut. Appart from the length L2 and the number of roots M2 = aL2 



^^It may be interesting to understand the relation to the string sphtting and joining approach of [46]. 

our normaHzation p{u), q{u) ^ 1 whereas J p{u)du — AI ^ L. This impHes that the argument of the 
exponent scales as L. 

^^For some configurations with large filling fractions we might find some logarithmic singularities in the 
first integral. The safest way to compute it in this case is to introduce some twists to regularize the integrals 
as done in [28] and then analytically continue the results to the case of zero twist. 
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n = 1 




Figure 5: In the classical limit the Bethe roots distribute themselves along disjoint cuts in 
the complex plane. The resolvent (100) has square root cuts at the position of the Bethe 
roots in the continuum limit and defines a Riemann surface. The final result for the structure 
constants (101) is given by some integrals in this surface. The integral involving the quasi- 
momenta q{v) is over the A-cycles of the surface denoted in green in the figure. The integral 
involving the density p{u) is taken over the cuts, n = 1, 2 are the mode numbers in this real 
data example. 



it has one more integer parameter - the mode number n which fixes the 27rin ambiguity of 
logarithm of the BAE (for more details see [27]). 

This solution is dual to a string state which, at strong coupling is described by a classical 
string motion, the rigid circular string. For concreteness let us briefiy describe this solution. 
The string is restricted to R x C AdS^ x S^, where R is the global AdS^ time t and 
the C can be describe by four embedding coordinates Xj satisfying ^^Xf = 1. The 
circular string classical solution is given by 



Wi V W2 

where a and r are the worldsheet coordinates in conformal gauge. The map to the spin chain 
state is established by L2 = ^J~\J and M2 = aV^jT". The others parameters of the solution 
are not all independent due to the equations of motion and the Virasoro constraints. The 
level matching condition gives m = an whereas the other constraints are more complicated 
and in the large J7 limit they give 

^ Q;(l — a)?T,^ ^ Q;(2a; — l)r;,^ ^ (a — l)(2a — l)n'^ , 
^"^^^ 2J ' "^^^^^ 2J . nj2^J+^- (104) 
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The classical string energy is K. 

Notice that J = L2I can be very large even for small couplings A ~ 0. This implies 
that the large J expansion of the classical string energy might capture some information 
about the weak coupling expansion. Indeed, expanding the anomalous dimensions in the 
scaling limit eq.(99) the one- loop spin chain gives 



72 = A- 



am 



2L, 



+ 



A 

-^2 



(105) 



which is precisely A — L^- This is a particular example of the FT limit alluded above. 

We hope that, similarly, the scaling limit of the structure constants eq.(lOl) can be 
obtained from the classical string theory in a similar limit. 

For illustration let us use the general expression (101) to evaluate C°23 in the small filling 
fraction a expansion. For simplicity we will focus on the case when M2 = iV2, i.e. the 
circular string is a highest weight, and expand for small a = M2/L2. It is useful to introduce 
/3 = (Li — N^)/L2 and r = nnp. We find (see appendix C for more details) 



l'-'123 I 



123 I ^^P 



L2T + O 



-^2 



where 



a log 



/ sinr^ 








V ^ J 





3/3^ 



+ 



(/3 



/3 



sin^ r 



(106) 



(107) 



5 Conclusions, speculations and open problems 

Computing three-point functions m. M = ^ super Yang-Mills is an ambitious goal with far 
reaching consequences. Together with the two-point functions, the three-point functions 
are enough to reconstruct any higher-point function in this conformal field theory. It is 
fascinating to note that integrability techniques, so useful in the spectrum problem, can be 
used to tackle the computation of structure constants at weak coupling. 

Basically, for the spectrum problem, we needed to know the spectrum of integrable Hamil- 
tonians. For the structure constants we require the eigenvectors as well, i.e. the precise 
structure of the spin chain wave functions. Using integrability we computed all tree-level 
non-extremal three-point functions of SU (2) sector operators, see table 2. The method used 
involved cutting and gluing single trace operators in a very stringy operation, see figure 1. 
It would be fascinating if we could interpret our results in a string field theoretic framework 
as advocated in [3]. The BMN limit might be a good starting point [29]. 

There are many interesting directions to pursue. 

We want to generalize our results to the full field strength multiplet of A/" = 4 SYM. After 
all, to glue three-point functions together into higher-point correlators we will need the most 
general cases. To complete this task we need to develop new technology for computing scalar 
products of Bethe states for higher rank groups such as PSU{2, 2|4). We will present a more 
detailed study of the Nested cases elsewhere [30]. Let us simply anticipate a few results. 
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As a first step we generalized the formula for the normalization of Bethe states to generic 
Lie (super) algebras r. We found a remarkably simple formula generalizing (66) which 
depends neatly on the Cartan matrix Ma^ and on the Dynkin labels Va of the corresponding 
algebra and representation. In that case the relevant phases are given by^^ 



= - log 




,=1 ^ - uP + IM,,, 



(108) 



where Ua'^ are the several Bethe roots of the corresponding Nested Bethe ansatz and there are 
Kj roots of each type where j = 1, . . . , r. For example, for the fundamental representation we 
have Vj = 6j^i. The normalization of the coordinate Bethe wave functions, with a convention 
analogue to (14) simply reads [30] 

r 

- ,,.a,U;;I,n-.,., deta,fe. (109) 

where we combined the indices a and j into a single index /. 

For the other rank one sectors of TV = 4, i.e. for the SL{2) and the SU{1, 1) sectors, we 
went a bit further. For these cases we generalized all possible scalar products (and not only 
the norm); this is discussed in section A. 2. With these results we can trivially convert the 
results of the main text to their analogue in the other rank one sectors.'^*' 

Another possible direction to generalize to more generic operators might be to consider 
the superprotected correlators of [33] as a vacuum and study excitations around these. 

One should also move to higher loops. There are two type of loop corrections. One type 
are loop corrections to the contraction of the three operators. At low loop orders, these are 
local corrections which dress the nonlocal overlap of the three wave functions which is studied 
in this paper. At one loop, and in the 50(6) sector, they are captured by Hamiltonian density 
insertions at the spin chain breaking points [3, 4, 5]. The insertion measures the energy cost 
of splitting the three operators [3]. 

The other type of corrections are loop corrections to the wave functions of the three 
operators. These comprise nonlocal corrections - due to the dependence of the scattering 
matrices (9) on the coupling - and also local corrections due to the insertion of the so-called 
fudge factors, or contact terms [34]. Due to the degeneracy of the spectrum at tree level, 
the computation of these corrections to the structure constants at I loops involves the (/ + 1) 
loops wave functions. 

In appendix E we list some basic examples of one-loop structure constants computed by 
taking into account these two types of corrections. 



^^To our knowledge the only previous results for scalar products of Bethe states in Nested Bethe ansatz 
systems are [31] for the norm of Bethe eigenstates in the SU{3) spin chain and [32] for the norm of Bethe 
eigenstates in the Hubbard model. We could not find any result for more generic inner products in the 
literature. 

^•^For example, we can go from the SU{2) to the SL{2) sector by formally replacing L —L in the several 
formulae (as usual at weak coupling). 
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Both loop corrections mentioned above could of course be incorporated in our treatment 
in a systematic way. In the meantime, given the experience with the spectrum problem, 
we could also try to make some educated guess from the analytic expressions for the one- 
loop structure constants. In the SU{2) spectrum problem, up to four loops, all one needed 
to modify in the asymptotic Bethe equations was the expression for the energy and the 
momentum of the excitations, e.g. 

e^^ = ^ + ^/^ ^ x{u + z/2) 

u-i/2 x{u-i/2) ^ ' 

where 



^x(«) = «+^«^-A (111) 

is the so-called Zhukoswky variable, introduced in [35]. If we look at our final expressions, 
e.g. at (97), we see that the building blocks are the same building blocks arising in the 
spectrum problem. It is therefore very tempting to try the simple replacement (110) and 
check whether our expressions could hold to higher loop order. 

By comparing the results obtained in this way with the brute force results of table 3 we 
conclude that such simple replacement does not work. 

On the other hand, given our discussion of the two type of loop corrections this is per- 
haps not surprising. Our treatment ignores contact terms and loop corrections to the Wick 
contractions. What we expect is that the replacement (110) captures the correct result for 
large (and maybe dilute) spin chains where the correction to the S'-matrix is the dominant 
effect. In contradistinction, the examples in the table 3 involved small operators. It could 
be that something a bit more sophisticated than the naive replacement (110) works even for 
small operators. We are currently investigating both possibilities. 

We should also try to explore the classical limit of the structure constants to make contact 
with strong coupling computations using the dual string sigma model. In the Frolov-Tseytlin 
(FT) limit 

A,J^oo with \' = ^<^l (112) 

the strong coupling string energy admits an expansion in A' which resembles a weak coupling 
expansion [38]. It turns out that the first two loop orders of the weak coupling anomalous 
dimensions are reproduced by this expansion. The match of weak and strong coupling 
results in the FT limit is most elegantly done through the algebraic curve method developed 
by KMMZ [39] . In this work it was explained that any string classical motion can be mapped 
to a Riemann surface defined by some quasimomenta q{x). The same is true at weak coupling: 
each solution to Bethe equations in the classical limit is described by some quasimomenta 
q{u). The main difference between the weak and the strong coupling quasi-momenta is in 
their pole singularities: in properly chosen variables the strong coupling quasimomenta has 
two poles at a; = ±A' while the weak coupling quasi-momenta has a single pole at m = 0, see 
(102). In the FT limit the strong coupling poles go to zero and the weak and strong coupling 



^^It is now understood that there is no deep reason for such match given the order or limits involved. Still 
it is a lucky coincidence that it works and we might also be fortunate to find it for the three-point functions. 
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Figure 6: AdS classical string solutions describing a three-point function in the CFT. (a) 
A correlation function between two heavy operators and one light operator is depicted. For 
examples of such computations at strong coupling see [36, 37]. (b) A decay of a heavy state 
into two light operators is depicted. Our prediction for the structure constant in this case is 
given by (97). It would be great if a simple upgrade of this expression a la KMMZ would 
produce the correct string theory value. In particular, it would be very interesting to study 
this process directly at strong coupling and check whether it matches (97) in the FT limit. 



curves coincide. This establishes a precise match of the spectrum in the FT limit for any 
classical state [39]. 

Optimistically, the same might be true for three-point functions. It would be great if the 
strong coupling result in the FT limit reproduced our classical result (101) which holds for 
any classical state. It might be simpler, at strong coupling, to consider more restricted string 
states. In (106) we presented the outcome of our result for the structure constant involving 
the operator dual to the rigid circular string and two BPS states. Correlation functions of 
one heavy operator and two light operators were not yet studied at strong coupling, see figure 
6b. Needless to say, it would be extremely interesting to do so. 

On the other hand, at strong coupling, there has been a lot of recent activity studying 
the correlation function of two heavy states and one light state [36, 37], see figure 6a. These 
works followed related earlier works on two-point functions [40] . It would be very interesting 
to take the continuum limit of our expressions in that case and compare with the strong 
coupling results (in the FT limit). We will report on a much more detailed study of the 
classical limit of our results elsewhere [23]. 
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A Scalar products 

In this appendix, we write several useful formulas for the scalar product Sn{{v}, {u}) for 
SU (2) spin chains in both the algebraic and coordinate bases. These complement the for- 
mulas in section 3. We also give explicit expressions for the scalar products for SL{2) and 
SU spin chains. In order to simplify the expressions, it is convenient to introduce the 
following functions 

h{u) = —— = —— , t{u) = —— = — — — , 113 
g{u) I f{u) u{u + i) 

and make use of the shorthand notation introduced in (43). 



A.l Scalar products in SU{2) 

General scalar product as a sum over partitions 

If all rapidities involved in the scalar product are arbitrary complex numbers, the scalar 
product can be written as a sum over all possible partitions of the two sets of rapidities. 
Explicitly [17, 41] 

SnUv}, {u}) = ^i"^^"*^ gi^^^"^ {-^)^"^^^'d'^a!'a?<rh''^K'^h'''^h?^ det det t^" , 

cxVJ a = {u\ 
7 U 7 = {v} 

(114) 

where we are using the notation introduced in (52) and (— 1)"^" is defined as a sign of the 
permutation of the ordered set {v} which gives aVJ a. Note that the sum runs over all 
partitions a U a = {u} and 7 U 7 = {v}, such that the number of elements in a and the 
number of elements in 7 are equal. We can use the conversion factor in (42) to obtain the 
corresponding expression in the coordinate base as: 

where is defined in (64). Recall that C{u*) = - [B{u)]'' so that 

s^^{{v*},{u}) = {-ir^\{v}\{u}r. (116) 
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New recursion relation in the coordinate basis 



Alternatively, we can define the general scalar product recursively, as in (73). Using the 
conversion factor (42), the new recursion relation in the coordinate basis reads 

Sn {{^U ■ ■ -^Vn}, {Ul, . . .,Un}) = ■■■,Vn,.. • , ^Tv}, ^2, • • ■ , MTv}) 

n 

- S^N-l{{uu Vi,...Vn,...,Vm,... Vn}, ^2, • • • , Un}) , (117) 

with the coefficients being in this case 

N N 



N 



N 



,co i^" 



v„.] - — — . . . 11 /(^i - ^^i) 11 '5(w„ 



d{ui)a{vr^ 



j>n 



9{ui + i)9{vn - i) 



N 



g{ui - V. 

j<m 



Y\_SiVn,Vj) Y\ S{Vj,Vm) + W S {Vm,V j)W_S {Vj ,Vn 



j>n 



j>m 



j<n 



d{ui) g{ui + I) g{vn - |) g{v„, - |) H /(^i - ^i) 

m 

a{ui) g{ui - ^) g{ui - Vn) g{ui - Vm) H f{vj-ui) 

n^n,m 

where we used the usual SU{2) S-matrix S{u,v) = f{u — v)/ f{y — u). 



Usual recursion relation 

Finally, for completeness, we write the recursion relation for the general scalar product known 
in the literature [17] 



StHx),d{x)] 



N N 

Yg{ui - Vn) Wg{ui 



N 



W g{Vk - Vn)k^n 



(118) 



n=l 



X {a'ivn)d{ui)S'^_i [a{x)h{ui — x), d{x)h{x — Vn)] 

— a{ui)d{vn)S'^_i [a{x)h{vn — x), d{x)h{x — ui)] } , 

where the sets of magnons entering S'^_i are obtained from the set entering S*^ by omitting 
Ul and Vn- The difference between both sides of (118) lies in the action of the operators A{u) 
and Vi^u) on the vacuum. For S^, this action is simply given by a{u) and d{u), see (32). 
For the first S'^_^ in (118), this action is given by a{u)h{ui — u) and d{u)h{u — Vn), while 
for the second S'^_^ it is a{u)h{vn — u) and d{u)h{u — Ui). That is, this recursion relation 
uses somehow generalized objects at each recursion stage, see also discussion at the end of 
section 3.3. 



39 



Scalar product with a Bethe eigenstate 

If the set of rapidities {u} satisfies the Bethe equations (15), while {v} are arbitrary complex 
numbers, the scalar product simplifies to [41] 



S^Mv}, {u}) = g^:^^''^ det n{u„ v,) , (119) 

where 

n{u„Vk) = a{vk)t{u, - v,)h^^'^''^ - {-lfd{vk)t{vk - u,)h^^'-^^ . 

In order to derive (119), one starts with (114) and makes use of the fact that the Bethe 
equations (15) for {u} can be written as 

We refer the reader to [41] for the details of the derivation. 
A.2 Scalar products in SL{2) and SU 

We now give several formulas for the scalar products in the other rank-1 subsectors of A/" = 4 
SYM. All these formulas can be checked against the brute force computation in the coordinate 
basis using the corresponding wave function. The latter is constructed in exactly the same 
way as in the SU{2) case, see (14), except that now the plave wave coefficients are obtained 
from products of the following S-matrices:'^^ 

Of ^ Ui) — Ua — t ^ 

^SL{2)[Ub, Ua) = , '-'5C/(1|1) = --■-• 

Ub-Ua + l 

SL{2) 

The general inner product for SL{2) spin chains obeys the same recursion relation as that 
for SU{2) (73), but with the following expressions for the building blocks 

a{u) = iu-]A , d{u) = iu+^- 



f(u) = , g{u) = - , h{u} = —— , t{u) = — — — . 

u u t u[u + l) 

■^^Recall that the SL{2) spin chain is non-compact, meaning that there is no bound on the number of 
excitations that we can have at a given site in the spin chain. That is, a generic SL{2) state is of the form 

|*)sL(2) = X! ^sL(2){ni,..-,nN)\ni,...,nN) . 

It is important to take into account these new limits of summation when computing the scalar product. 
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However, the factor that relates the general inner product in the algebraic and coordinate 
bases is different than that for SU{2). Namely 

Recall that the SL{2) Bethe equations are 

e^^^=l where ^ fn_^\' ^u^-u,^ _ 

\uj-i/2j l^.Uj-Uk-i 

Then, if the set of rapidities {u} satisfies these Bethe equations, the scalar product simplifies 
to 

Stm, {u}) = gi^^^"^ d^-^ det V,) , (122) 

where the matrix Vt is 

Finally, the norm in the coordinate basis is given by our conjectured formula (109) 

Ko = ,^^1, det d,<P, , (123) 

with (pj defined in (121). 



SU{1\1) 

The general inner product for SU{1\1) spin chains can be computed using the sum over 
partitions formula of SU{2) (114), but with the following expressions for the building blocks 

a{u) = + , d{u) = (^u-^ 

1 i 1 1 

/(m) = -, g{u) = -, h{u) = -, t{u) = — 

U U I u 

In this case, the relation between the scalar product in the algebraic and coordinate bases is 
exactly the same as that for SU{2) (115). The SU{l\i) Bethe equations are 

e"f>^ = l where e''^^ = Z ^^+V^ y .^24) 

Then, if the set of rapidities {u} satisfies these Bethe equations, the scalar product simplifies 
to 

Stm, {u}) = gi^^^''^ rfW det niu„ v,) , (125) 

where the matrix Q is now given by 

n{uj,Vk) = (-i)^ [a{vk)t{uj - Vk) + d{vk)t{vk - Uj)] . 

Finally, the norm in the coordinate basis is again given by formula (109) with 0j defined in 
(124). 
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B The BMN limit 



In this section we consider a simple limit when the lengths of operators La scale to infinity 
whereas the numbers of magnons Ma are fixed. We also assume that the momenta of the 
magnons are small (or equivalantly the Bethe roots are large). That is, similarly to the ther- 
modynamical limit we have Uk — La- These states describe small quasi-classical fiuctuations 
about the BMN point-like string. For related papers, see [7, 42]. 

In this limit in the Bethe ansatz equations one can replace the S matrix factor by 1 which 
leads to the trivial quantization condition: 



L 



Uk 



(126) 



where is an integer. In general one should assume that all n^'s are different. Otherwise 
there is a degeneracy which is lifted at the next order only."^'^ 

As it was shown in the main text, in order to construct all the structure constants one 
needs only tree key structures denoted by They can be easily expanded in the near 

BMN limit^^ 



^(/|M) 

^({«}) 
^(/|M,W) 



n(^ 

k=l 



L-2M 
N -M 



in — 



En 

o- k=l 



a 



a 



(128) 
(129) 
(130) 



Having these quantities at hand we simply combine them into structure constants using the 
expressions in table 2. For example. 



T-rMi v^^ T-rM2 
•'■•'■A:=l llfc=l 



M2 



'-^123 — 



. N-, -M 



I )\N2-M2 )\N^) 



L1L2-L3 



E E 



«12 "12 

e "fc — e'"'k 



V. 



e"k — e"k 



k=l 



(131) 



Namely, for the M„ roots with rik = n to the next order in are [43] 

^ - 1,...,M„ 



^— (l + iziV2L 
2Tm\ 



(127) 



where zi are the zeros of the Hermite polynomial, ilM^izk) = 0. In what follows we assume that the roots 
are large and well separated between each other. 

■^^It is also very simple to derive them directly from their definition and the form of the wave function 
(14). Recall that up to trivial factors SS^s^ are related to the norm of eigenstates, inner product with 
vacuum descendants and general inner product respectively. When the S'-niatrix is replaced by 1 all these 
quantities are very simple to compute directly without any fancy integrability machinery. 
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C Circular string details 



The density of roots in the scahng hmit and the resolvent are known exphcitly [39] 
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p{u) = ^^{27mu - Le+^<t'){27inu - Le-''t>) , (132) 

G(u) = I -^^dv = —J(27mu - Le+^<^)(2™ - Le"*"^) + — - nn (133) 
J u — V 2u 2u 

where the fiUing fraction M2/L2 = a is related to by a = sin^ |. Substituting these 
expressions into eq.(lOl) we get 

L2 SS2 2 \ ar^ J \ 6p^ 4 J 

3 (r^ (1 - 15(/3 - l)2/32(cos(2r) + 2) csc4(r)) (/3 - l)(2/3 - l)r^ cot(r) csc2(r) 1 

4 (/3 - 1)^(2/3 - l)r^(ll cos(r) + cos(3r)) 
~24 ~" 12sin^(r);g3 



(63(;9 - 1)3/33(26 cos(2r) + cos(4r) + 33) csc6(r) + 8) 

4536/36 

^r^ (120(/3 - l)/3(5(/3 - l)/3 + l)(cos(2r) + 2) csc4(r) - 8) 

1440/34 

where (3 = I12/L2 and r = 7rn/3. 



D Mathematica codes 

In this appendix we provide some Mathematica codes for computing the structure constants 
described in the main text. We implement both the brute force computation (85) as well as 
the final analytic results of table 2. Some examples are presented at the end of this appendix. 



Structure constants by brute force (85) 

Of f [Det : :matsq] ; Le = Length; 

f[u_] = 1 + I/u; g[u_] = I/u; n[0] = 0; 

S[x_, yj := (x - y + I)/(x - y - I) 

Wave [l_List] := Block [{p = Permutations [Range [Le [1] ]] , i, j}. Sum [A [p [ [i] ] ] 
Product [((l[[p[[i, j]]]] + I/2)/(l[[p[[i, j]]]] - I/2))-n[j], {j , 1, Le [1] }] , 
{i, 1, Le[p]>] //. {A[{a___, b_, c_, d }] :> S[l[[b]], l[[c]]] A[{a, c, b, d>] 

^^The algebraic curve for the general SU (2) circular string at strong coupling was similarly written down 
in [44]. 
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/; b > c} /. {A[a___] :> 1 /; a == Range [Le [a] ] }] ; 

normbf[L_, l_List] := Sum [ (Wave [1] /. Complex [a_, b_] -> a - I b) Wave [1] , 
Evaluate [Sequence M Table [{n[j], n[j - 1] + 1, L}, {j , Le[l]}]]] 

dphi [L_ , l_List] : = DetQTable [-If [i==j ,L/ (1 [ [i] ] "2+1/4) -Sum [2/ ( (1 [ [i] ] -1 [ [k] ] ) "2+1) , 
{k,Le[l]}] ,0]-2/(l+(l[[i]]-l[[j]])-2),{i,Le[l]},{j,Le[l]}] /. Det [{}] -> 1 

pref actor [l_List] := Product [l/(g[l[[j]] + I/2]g[l[[j]] - 1/2]), { j , 1, Le [1] }] * 
Product [f[l[[j]] - l[[i]]]/f [Conjugate [l[[j]] - l[[i]]]], {i, 1, Le[l]}, 
{j, i + 1, Le[l]>] 

normdet[L_, l_List] := pref actor [1] dphi[L, 1] 

C123[L1_, Nl_, L2_, N2_, L3_, N3_, ll_List, 12_List, 13_List] := 
Block [{i, j, psis, norms, limits}, 

psis = (Wave [11] /. n[j_] :> LI - N3 + j - N2 /; j > N2) 

(Wave[12] /. n[j_] -> L2 + 1 - n[N2 - j + 1] ) (Wave [13] /. n[j_] -> j); 

norms = normdet[Ll, 11] normdet[L2, 12] normdet[L3, 13]; 

limits = Sequence M Table[{n[j], n[j - 1] + 1, LI - N3>, { j , N2}] ; 

Sqrt [LI L2 L3/norms] If [limits === Sequence [] , psis. Sum [psis, limits]]] 

Note that the norms can be computed by brute force using normbf, as in (86), or using 
normdet, which is just the implementation of (66). We choose to use the latter as it is 
computationally much more efficient. 

Structure constants from table 2 

Of f [Det : :matsq] ; Le = Length; 

f[u_] = 1 + I/u; g[u_] = I/u; h[u_] = f[u]/g[u]; t [uj = g[u]"2/f[u]; 
f[ll_List, 12_List] := Product [f [11 [ [j 1] ] - 12 [[j2]]], {jl, Le[ll]}, {j2, Le [12] }] 
h[ll_List, 12_List] := Product [h [11 [ [j 1] ] - 12 [[j2]]], {jl, Le[ll]}, {j2, Le [12] }] 
fs[ll_List] := Product[f [ll[[jl]] - ll[[j2]]], {jl, Le[ll]}, {j2, jl + 1, Le[ll]}] 
gs[ll_List] := Product [g [11 [[jl]] - ll[[j2]]], {jl, Le[ll]}, {j2, jl + 1, Le[ll]}] 
fb[ll_List] := Product[f [ll[[jl]] - ll[[j2]]], {jl, Le[ll]}, {j2, jl - 1}] 
gb[ll_List] := Product [g [11 [[jl]] - ll[[j2]]], {jl, Le[ll]}, {j2, jl - 1>] 
gp[l_List] := Times 00 g[l + 1/2]; gm[l_List] := Times @@ g[l - 1/2] 
a[l_List] := Times m ((1 + I/2)/(l - 1/2)) 

sign[a_List, ab_List, v_List] := Signature [Join [a, ab]] Signature [v] 

Dvd[ls_List] := ({Complement [Is , #1], #1} & ) /@ Subsets [Is, {0, Le [Is] }] ; 

dphi [L_ , l_List] : =Det®Table [-If [i==j , L/ (1 [ [i] ] "2+1/4) -Sum [2/ ( (1 [ [i] ] -1 [ [k] ] ) "2+1) , 

{k,Le[l]}] ,0]-2/(l+(l[[i]]-l[[j]])"2),{i,Le[l]},{j,Le[l]}] /. Det [{>] -> 1 
dett [ll_List , 12_List] : =Det@Table [t [11 [ [i] ] -12 [ [ j ] ] ] , {i , Le [11] } , 

{j,Le[12]}] /. Det[{}] -> 1 
A[L_, ls_List] := Block [{dv = Dvd[ls] ,al,alb}, Sum[al = dv[[i, 1]]; 

alb = dv[[i, 2]]; (-l)"Le[al] f[al, alb] /a[al] "L, {i, Le[dv]}]]; 
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B [L_ , N_ , l_List] : =gni [1] f s [1] Sqrt [f b [1] dphi [L , 1] Binomial [L-2Le [1] , N-Le [1] ] 
/ (f b [Conjugate [1] ] L) ] /Sqrt [gp [1] gm [1] ] 

T [n_ , u_List , v_List] : =Block [{dv=Dvd [v] , du=Dvd [u] , al , alb , be , beb} , gs [u] gb [v] Sum [ 
al=dv [ [i , 1] ] ; alb=dv [[1,2]]; be=du [ [ j , 1] ] ; beb=du [ [ j , 2] ] ; If [Le [al] ==Le [be] && 
Le [alb] ==Le [beb] , sign [al , alb, v] sign [be , beb, u] dett [be,al] dett [alb, beb] a [beb] "n 
a [al] "n h [be , al] h [alb , beb] h [be , beb] h [alb , al] , 0] , {i , Le [dv] } , { j , Le [du] }] ] 

Cooo [Ll_ , Nl_ , L2_ , N2_ , L3_ , N3_] : = (Binomial [L1-N1+N2 , N2] / 

(B[L1,N1,{>]B[L2,N2,{>]B[L3,N3,{>])) 
Cxoo [Ll_ , Nl_ , L2_ , N2_ , L3_ , N3_ , ll_List] : = (Binomial [Ll-Nl+N2-Le [11] , N2] 

A [N1-N2 , 11] / (B [LI , Nl , 11] B [L2 , N2 , {}] B [L3 , N3 , {>] ) ) 
Coxo [Ll_ , Nl_ , L2_ , N2_ , L3_ , N3_ , 12_List] : = (Binomial [Ll-Nl+N2-Le [12] , N2-Le [12] ] 

A [L1-N1+N2 , 12] / (B [LI , Nl , {}] B [L2 , N2 , 12] B [L3 , N3 , {}] ) ) 
Coox [Ll_ , Nl_ , L2_ , N2_ , L3_ , N3_ , 13_List] : = (Binomial [L1-N1+N2 , N2] 

A [N1-N2 , 13] / (B [LI , Nl , {>] B [L2 , N2 , {}] B [L3 , N3 , 13] ) ) 
Coxx [Ll_ , Nl_ , L2_ , N2_ , L3_ , N3_ , 12_List , 13_List] : = (Binomial [Ll-Nl+N2-Le [12] , N2-Le [12] ] 

A [L1-N1+N2 , 12] A [N1-N2 , 13] / (B [Li , Nl , {}] B [L2 , N2 , 12] B [L3 , N3 , 13] ) ) 
Cxox [Ll_ , Nl_ , L2_ , N2_ , L3_ , N3_ , ll_List , 13_List] : = (Binomial [Ll-Nl+N2-Le [11] , N2] 

A [N1-N2 , 11] A [N1-N2 , 13] / (B [LI , Nl , 11] B [L2 , N2 , {}] B [L3 , N3 , 13] ) ) 
Cxxo [Ll_,Nl_,L2_,N2_,L3_,N3_,ll_List,12_List] : = (dv=Dvd [11] ; 

Sum [al=dv [ [i , 1] ] ; alb=dv [ [i , 2] ] ; If [Le [al] ==N1-N2 , a [al] "LI f [alb , al] A [N1-N2 , al] 

T [L1-N1+N2 , alb , 12] , 0] , {i , Le [dv] }] / (B [LI , Nl , 11] B [L2 , N2 , 12] B [L3 , N3 , {>] ) ) 
Cxxx [Ll_ , Nl_ , L2_ , N2_ , L3_ , N3_ , ll_List , 12_List , 13_List] : = (A [N1-N2 , 13] B [L3 , N3 , {}] 

Cxxo[Ll,Nl,L2,N2,L3,N3,ll,12]/B[L3, N3, 13]) 

Examples 

Let us show how to use the code above to compute the most general case C^** in table 2 for 
some three-point function configuration. In this example, the sets of rapidities us, vs, ws 
satisfy the Bethe equations of operators Oi, O2 and O3, respectively, with the lengths and 
number of excitations indicated in each case. Evaluating 

Ll=9; Nl=4; L2=9; N2=2; L3=4; N3=2; 

us = {0.414080361016 - 0.993048580811 I, 0.409292874229, 

0.414080361016 + 0.993048580811 I, -0.131911999475}; 
vs = {-0.207106781187, 0.207106781187}; 
ws = {-0.288675134595, 0.288675134595}; 
C123 [Ll,Nl,L2,N2,L3,N3,us,vs,ws] 
Cxxx [L 1 , Nl , L2 , N2 , L3 , N3 , us , vs , ws] 

we obtain a perfect match between the brute force computation (85) and the analytic pre- 
diction of table 2. For these values of the Bethe roots we find 

C"' = -0.473631 + 0.079146 i . 
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We can also use the codes above to get analytic results. For example, to reproduce C^'^ 
for the case of two magnons N2 = 2, see equation (98), we would simply run 

ClearAll [LI , Nl , L2 , N2 , L3 , N3] 

f si=FullSimplif y [# , {0<p<\ [Pi] /2}] & ; 

fsi2=#//.{Exp[a_]-Exp[b_] :>Exp[(a + b)/2]2 Sinh[(a-b)/2//fsi] , 

Exp[a_]-l:>Exp[a/2]2 Sinh [a/2//f si] ,Exp[a_] : >Exp [a//f si] >& ; 
f siSCoxo [LI , N3+2 , L2 , 2 , L3 , N3 , {u , -u}] / . u->Cot [p/2] /2//TrigToExp//Factor//f si//f si2 



E Data 

In this appendix we will produce one-loop data for the future. Two things need to be done: 
we need to contract the external wave functions of the three operators using Hamiltonian 
insertions [3, 4, 5] and we need to take into account the C(A) correction to the wave function, 
which is a two- loop effect. The latter effect was not taken into account in previous works 
[3, 4, 5] but its importance was mentioned in the conclusions of [3]. 

For simplicity we shall consider two kind of operators only: BPS operators and operators 
with two impurities with opposite momentum 

/I 7Tn XL 2Tm ,^,,9x\ , s. 

" = ±(2"'r3T + s5r3T=-rTi+o(^^)j (135) 

and energy 

A 9 riTT }? \ ( ^ 27m \ . irn ^^/^^x 

^ = ^ r^T - 4^ 131 (1 + i + 2 r^l) 131 + 0(A») . (136) 

The integer n = l,...,L/2 — 1 is called the mode number. These states are denoted by 
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{n,L) and they diagonalize the two-loop Hamiltonian of [45]. We take'^^ {g'^ = A/IGtt^) 

(l,4) = Tr[Z,X]' , 

(1, 5) = Tr {Z^X^) - Tr {Z^XZX) , 



(1,6) = Tr (Z^X^) 
(2, 6) = Tr {Z^X^) 



^ "^(1 + g') Tr (Z'XZX) + ^ v^+(3 v^^ (Z^XZ^X) 



2 

3 + 



0^ 1 

(1, 7) = Tr [Z^X^] - y Tr [Z^XZX) ~ ^{2 - Tr [Z^XZ^X] , 

(2, 7) = Tr (Z^X^) - (^2 + ^1^^ Tr (Z^XZX) + (^1 + ^^ Tr (Z^XZ^X) , 



Tr (Z^X^) + 



7 



TT \ 



sec - sin — - 0.489997 g 



7r\ / Stt 
sec ( — sm 



77 \U 



UJ 



0.206966 g^ 



Tr (Z^XZX) 



Tr (Z^XZ^X) 



1 /TT 



2 sec (^y ) - 0.283031 g' 



(2, 8) = Tr (Z^X^) 



cos 



+ 



TT 

sm I — ) CSC 
14 

371 

14 



IT 

7 

371 

14 



Tr (Z^XZ^X) , 
37r\ 



CSC ( ^ 1 + 1.117057^2 



1.023185 g^ 



Tr (Z^XZX) 
Tr {Z^XZ^X) 



CSC — + 0.093872 g' 



Tr (Z^XZ^X) 



(3, 8) = Tr [Z^X^] 



vr \ /37r 
CSC I — sin 



14/ V 14 



+ 3.392947/ 



Tr {Z^XZX) 



+ 



/ TT \ 



TT 



CSC j COS (^y ) + 8.769849 g 
'1 / vr \ 



Tr {Z^XZ'^X) 



2 CSC (^-j +5.376903 



Tr {Z'^XZ^X) . 



Note that for convenience, we gave the numerical values of exact expressions in the states 
shown above. For example, in the last term of 03,8, we have 



5.376903 



9-8cos(f) +8sin(ff) 



44 - 54 cos (f) +78 sin (^) -20sin(ff) 



■^^In this section we shall not keep track of the overall phase of the structure constant but only its absolute 
value. 
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O2 


O3 


r 




O2 


O3 


r 


(1,4) 


BPS 


BPS 


-6 


(3,8) 


BPS 


BPS 


-9.561263 


BPS 


BPS 


(1,4) 


-6 


[2,4] 


(3,8) 


[0,4] 


-23.396074 


(1,5) 


BPS 


BPS 


-4 


[4,4] 


(3,8) 


[2,8] 


-9.561263 


BPS 


BPS 


(1,5) 


-4 


[4,8] 


(3,8) 


[2,4] 


-9.561263 


(1,6) 


BPS 


BPS 


-11/2 + 13\/5/10 


[3,5] 


(3,8) 


[1,5] 


-8.188323 


[3,4] 


(1,6) 


[1,4] 


3a/5/10 - 7/2 


[4,6] 


(3,8) 


[2,6] 


-8.188323 


[4,5] 


(1,6) 


[2,5] 


3^/5/10 - 7/2 


[3,6] 


(3,8) 


[1,4] 


-5.330327 


BPS 


BPS 


(1,6) 


-11/2 + 13^5/10 


[4,7] 


(3,8) 


[2,5] 


-5.330327 


(2,6) 


BPS 


BPS 


-11/2-13^/5/10 


[4,5] 


(3,8) 


[2,7] 


-5.330327 


[3,4] 


(2,6) 


[1,4] 


-3^5/10-7/2 


[3,4] 


(3,8) 


[1,6] 


-5.330327 


[4,5] 


(2,6) 


[2,5] 


-3^5/10-7/2 


BPS 


BPS 


(3,8) 


-9.561263 


BPS 


BPS 


(2,6) 


-11/2-13^5/10 


BPS 


(1,4) 


(1,4) 


-12 


(1,7) 


BPS 


BPS 


-7/4 


[4,5] 


(1,5) 


(1,4) 


-10 


[3,4] 


(1,7) 


[1,5] 


-2 


(1,4) 


[0,5] 


(1,5) 


-10 


[3,5] 


(1,7) 


[1,4] 


-2 


[4,6] 


(1,6) 


(1,4) 


13^/5/10-23/2 


[4,5] 


(1,7) 


[2,6] 


-2 


(1,4) 


[0,6] 


(1,6) 


13\/5/10- 23/2 


[4,6] 


(1,7) 


[2,5] 


-2 


[4,7] 


(1,7) 


(1,4) 


-31/4 


BPS 


BPS 


(1,7) 


-7/4 


(1,4) 


[0,7] 


(1,7) 


-31/4 


(2,7) 


BPS 


BPS 


-27/4 


[4,8] 


(1,8) 


(1,4) 


-7.237755 


BPS 


BPS 


(2,7) 


-27/4 


(1,4) 


[0,8] 


(1,8) 


-7.237755 


(1,8) 


[0,8] 


[2,4] 


-1.237755 


(1,4) 


(1,8) 


[0,4] 


5.123727 


[2,4] 


(1,8) 


[0,4] 


-4.525727 


[4,6] 


(2,6) 


(1,4) 


-13^5/10-23/2 


[4,6] 


(1,8) 


[2,6] 


-1.513645 


(1,4) 


[0,6] 


(2,6) 


-13^5/10-23/2 


[3,5] 


(1,8) 


[1,5] 


-1.513645 


[4,7] 


(2,7) 


(1,4) 


-51/4 


[3,4] 


(1,8) 


[1,6] 


-1.455824 


(1,4) 


[0,7] 


(2,7) 


-51/4 


[3,6] 


(1,8) 


[1,4] 


-1.455824 


[4,8] 


(2,8) 


(1,4) 


-11.200983 


[4,5] 


(1,8) 


[2,7] 


-1.455824 


(1,4) 


[0,8] 


(2,8) 


-11.200983 


[4,7] 


(1,8) 


[2,5] 


-1.455824 


(1,4) 


(2,8) 


[0,4] 


-1.554117 


[4,4] 


(1,8) 


[2,8] 


-1.237755 


[4,8] 


(3,8) 


(1,4) 


-15.561263 


[4,8] 


(1,8) 


[2,4] 


-1.237755 


(1,4) 


[0,8] 


(3,8) 


-15.561263 


BPS 


BPS 


(1,8) 


-1.237755 


(1,4) 


(3,8) 


[0,4] 


-9.293748 


(2,8) 


[0,8] 


[2,4] 


-5.200983 


[4,4] 


(1,5) 


(1,5) 


-8 


[2,4] 


(2,8) 


[0,4] 


-10.078199 


(1,5) 


[0,6] 


(1,5) 


-8 


[4,7] 


(2,8) 


[2,5] 


-7.213849 


[4,5] 


(1,6) 


(1,5) 


3^5/10-15/2 


[4,5] 


(2,8) 


[2,7] 


-7.213849 


(1,5) 


[0,7] 


(1,6) 


13\/5/10- 19/2 


[3,6] 


(2,8) 


[1,4] 


-7.213849 


[4,6] 


(1,7) 


(1,5) 


-6 


[3,4] 


(2,8) 


[1,6] 


-7.213849 


(1,5) 


[0,8] 


(1,7) 


-23/4 


[4,8] 


(2,8) 


[2,4] 


-5.200983 


[4,7] 


(1,8) 


(1,5) 


-5.455824 


[4,4] 


(2,8) 


[2,8] 


-5.200983 


[4,5] 


(2,6) 


(1,5) 


-3v^/10- 15/2 


[3,5] 


(2,8) 


[1,5] 


-0.298032 


(1,5) 


[0,7] 


(2,6) 


-13^5/10- 19/2 


[4,6] 


(2,8) 


[2,6] 


-0.298032 


(1,5) 


[0,8] 


(2,7) 


-43/3 


BPS 


BPS 


(2,8) 


-5.200983 


[4,7] 


(2,8) 


(1,5) 


-11.213849 



Table 3: Data for the ratio between the one-loop and tree-level structure constants parametrized 
by r in A^cC'123 = c^^]. (l + g^f + 0(9^)) • The notation (n, L) indicates a two magnon operator with 
mode number n and length L. The notation [A^, L] indicates a vacuum descendant with spin 
flips and total length L. The vacuum and excitation choice for each of the three operators is given 
in (81). When we write BPS we can replace it by any [N, L] with < 4 and L < 8 (such that the 
three-point function exists); in this case the result is independent of and L. To make the table 
easier to read we colored the BPS operators in red and the non-BPS operators in blue. We also 
colored in green the ratios for configurations with /13 = 0. For these cases (only for these cases) 
the mixing of O2 with double trace operators should in principle be included [7, 8] (we did not take 
this effect into account). 



We can also write any protected state with excitations and length L, which we denote by 
[N,L], as 



For example: 

[2, 5] = 5 Tr (Z^X^) + 5 Tr (Z^XZX) , 

[3, 6] = 6 Tr {Z'^X'^) + 6 Tr {Z^XZX^) + 6 Tr (Z^X^ZX) + 2 Tr (ZXZXZX) . 

All the states written in this appendix consider the vacuum as being Z fields and the 
excitations as being the X fields. In other words, they are good states for the operator Oi. 
For the operators O2 and O3 we use the same states but with different scalars playing the 
role of vacuum and excitations, see (81). 

Having computed the two-loop eigenstates (n, L) and knowing how to write any BPS 
operators with a given length and number of excitations, we provide a list of ratios of one- 
loop structure constants to tree-level structure constants in table 3 involving these two kind 
of operators. 
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